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ABSTRACT 

In the three-dimensional parameter space defined by velocity dispersion (c), effective radius (i? e ), 
and effective surface brightness (I e ), early- type galaxies are observed to populate a two-dimensional 
fundamental plane (FP) with finite thickness. In Paper III of this series, we showed that the thickness 
of the FP is predominantly due to variations in the stellar mass surface density (£*) inside the effective 
radius R e . These variations represent differences in the dark matter fraction inside R e (or possibly 
differences in the initial mass function) from galaxy to galaxy. This means that galaxies do not wind 
up below the FP at lower surface brightness due to the passive fading of their stellar populations; 
they are structurally different. Here, we show that these variations in £* at fixed dynamical mass 
(M,iyn) are linked to differences in the galaxy stellar populations, and therefore to differences in their 
star formation histories. We demonstrate that the ensemble of stellar population and £* variations 
through the FP thickness can be explained by a model in which early-type galaxies at fixed Md yn have 
their star formation truncated at different times. The thickness of the FP can therefore be interpreted 
as a sequence of truncation times. Galaxies below the FP have earlier truncation times for a given 
Mdyn, resulting in lower £*, older ages, lower metallicities in both [Fe/H] and [Mg/H], and higher 
[Mg/Fe] . We show that this model is quantitatively consistent with simple expectations for chemical 
enrichment in galaxies. We also present fitting functions for luminosity- weighted age, [Fe/H], [Mg/H], 
and [Mg/Fe] as functions of the FP parameters a, R e , and I e . These provide a new tool for estimating 
the stellar population properties of quiescent early-type galaxies for which high-quality spectra are 
not available. 
Subject headings: galaxies: evolution, galaxies: abundances, galaxies: elliptical and lenticular 



1. INTRODUCTION 

In the cold dark matter (CDM) paradigm, galax- 
ies form hierarchically. However, more massive galax- 
ies are composed of older stellar populations than are 
less massive gala xies (so-ca ll ed "a rcheological down- 
sizing", e.g., Thomas " et al.l I2005D . Modern semi- 
analytic models (SAMs) for galaxy formation have shown 
that archeological down-sizing can be produced within 
the curren t ACDM framework of concordance cosmol- 
ogy (e.g.. iDe Lucia et all 120061: jCattaneo et ail 120081: 
iSomerville et al.l2008tlWang et all2008D . The key to rec- 
onciling archeological down-sizing with hierarchical as- 
sembly is to recognize that the first process refers to the 
formation of stars, while the second process refers to the 
mass build-up of a galaxy, and that these two processes 
can be decoupled. At the same time, galaxies are an 
orderly family and obey numerous scaling relations be- 
tween their stellar populations and structural properties. 
The tightness of the relations indicates that star forma- 
tion in galaxies and the mass assembly process, though 
decoupled, cannot be combined at random. 

This is particularly true of the quiescent, non-star- 
forming galaxies discussed in this series of papers. To 
zeroth-order, early-type galaxies form a one-dimensional 
(ID) sequence, with their star formation histories and 
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many of their structural properties showing strong 
trends with galaxy mass. However, they are not a 
purely ID family: early- type galaxi es populate the two- 
dimensional Fundamenta l Plane (jDjorgov ski fe Davisl 
119871: iDressler et~al~lll987l FP) in the three-dimensional 
parameter space of galaxy properties defined by central 
velocity dispersion (a), effective radius (R e ), and effec- 
tive surface brightness (I e = L/2ttRI, where L is galaxy 
luminosity). In fact, the FP has finite thickness, which 
adds a third dimension to the early-type galaxy struc- 
tural parameter space. Thus the ID nature of many 
galaxy scaling relations disguises a multi-dimensional 
galaxy family. 

While the high-mass FP appears nearly unchanged 
since z ~ 1 other than expected passive evolution in 
M+/L, low- mass galaxies are still "settling" onto the F P 
since z ~ 1 (Ivan der Wei et al.ll200llTreu et al1l2005allb1 : 
Ivan der Wei et all 120051 ). The local FP should then be 
expected to include at least some galaxies that are rela- 
tiv ely recent ar ri vals. This is supported by the result s 
of iForbes et all (fl99l and iTerlevich fc Forbes! (|2002l ). 
who find that residuals from the FP correlate with age 
such that galaxies offset to higher (lower) surface bright- 
ness have younger (older) ages than those occupying the 
FP, based on simple stellar population (SSP, i.e., single 
burst) models. The thickness of the FP may therefore be 
an age sequence, with the most recent arrivals lying at 
higher surface brightness than the midplane. This thick- 
ness reveals a (narrow) third dimension to the structural 
parameters of early-type galaxies. 

This paper is the fourth installment in a series of pa- 
pers that seeks to connect the star formation histories 
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of early-type galaxies to their global structural parame- 
ters. We explore these relations in an explicitly multi- 
dimensional parameter space, which allows us to identify 
the fundamental trends in a model-free fashion. 

It has been well-documented that the stellar popu- 
lation properties of early-type gal axies scale with their 
stella r mass (M+) o r with a (e.g.. iKuntschner fe Daviei 
19981: iTraeer et all l2000aJ: iTerlevich fc Forbesl 12002 



Proctor fc Sansoirir 2002t Thomas et al 



2005 



NelanetalJ l2005t iGaJlazzi et all 120051 : ISmith et al.1 
2007t iGraves et al.1 I2007T) . However just as the Faber- 
Jackson relation (jFaber fc Jackson! I1976T ) turns out to 
be a ID projection of the 2D Fundament al Plane, we 
showe d in Papers I and II of this series (|Graves et al.l 
l2009allbT ) that the above-mentioned familiar stellar 
population trends are in fact ID projections of an 
underlying 2D family of stellar population properties. 

The main dimension consists of correlations be- 
tween stellar po pul ation pr operties and a 
Proctor fc Sansoml l2002t iThomas et ail 
Nelan et al.l l2005f) . One of our main findings in 
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Papers I and II was to demonstrate that in fact the star 
formation histories of galaxies depend most strongly 
on the velocity dispersion a and not on stellar mass, 
dynamical mass, or luminosity. Thus the various 
"mass-dependent" relations are in reality er-dependent 
relations. This conclusion has been corroborated by the 
result s of Ivan der Wei et all ([200 9) and iRogers et al.1 
(2010). In addition to this first dimension of variation, 
numerous authors have reported a spread in stellar 
population ages and metallicties at fixed a, with an anti- 
correlation observed between age and metallicity (e.g. 



Worthev et alj|1995l: IColless et al.lll9ffi[J0rgensenlll999l; 
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2008a). In Papers I and II, we not only confirm the 
existence of this second dimension of variation, we also 
show that it is strongly correlated with other global 
galaxy properties. 

Paper I mapped the 2D family of stellar population 
properties onto a familiar ID early-type galaxy scal- 
ing relation: the color-magnitude relation (CMR). We 
showed that the various stellar population properties be- 
have differently with respect to a and to magnitude, such 
that age and [Mg/Fe] show stronger trends with a, while 
[Fc/H] and [Mg/H] show stronger trends as functions of 
magnitude. A 3D mapping of stellar population proper- 
ties onto a, color, and magnitude demonstrated that the 
CMR hides a 2D family of stellar population properties. 
The first dimension of this family maps most strongly 
onto a, while the second dimension correlates with mag- 
nitude residuals at fixed <70 

Paper II carried the analysis to the standard 3D space 
of global galaxy properties: the Fundamental Plane space 
defined by er, i? e , and I e . This set of parameters is 
our preferred space in which to explore stellar popula- 
tion trends because of two important points shown in 
Paper II: (1) that on the FP midplane, galaxy star for- 



5 In a subsequent analysis, Smith et al. (2009) find a similar 2D 
family of stellar population variations in a sample of galaxies from 
the Shapley Supercluster. However, they interpret these as due to 
aperture effects in the observations. This interpretation is incon- 
sistent with the results demonstrated in Paper II, which strongly 
support our original interpretation that the stellar population ef- 
fects are intrinsic and genuine. 



mation histories are independent of R e and therefore de- 
pend explicitly on a, not on total galaxy mass, and (2) 
that the second dimension of stellar population varia- 
tion correlates strongly with surface brightness residuals 
at fixed a and R e . We showed that sorting galaxies by 
surface brightness residuals, instead of L residuals (as 
in Paper I) , gives a cleaner differentiation of the various 
stellar population properties, and thus that it is A/ e , not 
AL, that drives the second dimension of variation. This 
means that the second stellar population dimension maps 
onto a cross-section through the FP. This result was new. 
It demonstrated not only that the thickness of the FP is 
real, but that it is an important parameter in the evo- 
lutionary history of galaxies. Furthermore, these trends 
exist for galaxies at fixed physical size i? e , which makes 
it impossible that the tren ds could be due to aperture 
effects alone as claimed by ISmith et al.l (|2009f ) . Taken 
together, Papers I and II show that the 2D family of 
galaxy star formation histories maps smoothly onto 3D 
FP space. This paper quantifies and expands on that 
work; 

In lGraves fc Faben (|2010l . hereafter Paper III) , we com- 
puted the variations in M+/L implied by the above 3D 
mapping in FP space (assuming a constant initial mass 
function) and compared these to the measured dynam- 
ical mass-to-light ratios (Md yn /L). We showed that, if 
the IMF is constant, M+/L variations are too small to ac- 
count for the observed variations in Md yn / L, both along 
the FP midplane (the "tilt" of the plane) and through 
the thickness of the plane. Hence, either the IMF is 
changing throughout FP space in a systematic way or the 
fractional contribution of dark matter relative to stars is 
changing. Focusing on the variations through the thick- 
ness of the FP (i.e., variations with I e at fixed a and fixed 
R e ), we considered the implications of IMF variations 
and showed that the most natural IMF scenarios predict 
abundance and line index variations that are qualitatively 
at variance with those observed. This led us to tentatively 
conclude that our estimates of stellar mass are valid and 
further that the observed variations in Mdyn/L imply 
real changes in stellar mass surface density (£*) through 
the plane; galaxies below the FP (at lower I e ) must have 
fewer stars within i? e , not just dimmer stars. This im- 
plies that galaxies do not evolve below the FP by passive 
fading — they are built that way. We went on to consider 
ways this might happen and proposed a toy model in 
which objects with the same a all start forming stars at 
roughly the same time but have different durations of star 
formation. In this picture, galaxies below the plane stop 
forming stars sooner (i.e., are "prematurely truncated") 
and end up with relatively less stellar mass, while galax- 
ies above the FP form stars longer, have higher overall 
efficiencies in converting gas into stars, and wind up with 
more stellar mass. We showed that this hypothesis is 
qualitatively consistent with all structural and line-index 
trends seen through the FP. 

This paper develops the theme of trends through the 
thickness of the FP still further. Here we make the re- 
lation between stellar populations and structural vari- 
ables even more explicit by showing in detail how the 2D 
family of stellar populatio n variations— essentia lly the 
"metallicity hyperplane" of Tra ger et al.l (|2000aj) — map 
onto structural slices through the plane. The original 
metallicity hyperplane involves just SSP age and [Fe/H]. 



Table 1 

Properties of Galaxy Bins and Stacked Spectra 



a bin 

(km s _1 ) 



A log I e bin 

(£© pc- 2 ) 



Median 
logo- 

(km s- 1 ) 



Median 

A log I e 

(L Q pc- 2 ) 



N* 



S/JVt 

(A- 1 ) 



H/3 

(A) 



(Fe> 
(A) 



Mg 6 

(A) 



[Om]A5007* 

(A) 



Ages 
(Gyr) 



[Fe/H]§ 
(dex) 



[Mg/H]§ 

(dex) 



[Mg/Fc]§ 
(dex) 



1.86- 


-2.00 


-0.25 .. 


.-0.15 


1.96 


-0.19 


139 


175 


1.91±0.05 


2.40±0.07 


3.23±0.06 


-0.15±0.03 


8.1±0.8 


-0.26±0.04 


-0.16±0.05 


0.10±0.03 






-0.15 .. 


.-0.05 


1.95 


-0.09 


278 


259 


2.11±0.03 


2.33±0.05 


3.32±0.04 


-0.25±0.02 


5.5±0.3 


-0.22±0.03 


-0.04±0.04 


0.18±0.03 






-0.05 .. 


. 0.05 


1.94 


0.00 


399 


332 


2.16±0.03 


2.37±0.04 


3.27±0.03 


-0.19±0.02 


5.0±0.2 


-0.18±0.02 


-0.03±0.03 


0.15±0.02 






0.05 .. 


. 0.15 


1.95 


0.10 


313 


318 


2.20±0.03 


2.36±0.04 


3.22±0.03 


-0.19±0.02 


4.6±0.3 


-0.17±0.02 


-0.03±0.03 


0.14±0.02 






0.15 .. 


. 0.25 


1.94 


0.18 


107 


201 


2.40±0.04 


2.37±0.06 


3.06±0.05 


-0.20±0.03 


3.1±0.2 


-0.10±0.03 


0.02±0.04 


0.12±0.03 


2.00- 


-2.09 


-0.25 .. 


.-0.15 


2.05 


-0.18 


216 


241 


1.96±0.04 


2.41±0.05 


3.42±0.04 


-0.25±0.02 


7.3±0.5 


-0.22±0.03 


-0.08±0.04 


0.14±0.03 






-0.15 .. 


.-0.05 


2.05 


-0.09 


556 


111 


2.09±0.02 


2.42±0.03 


3.50±0.02 


-0.28±0.01 


5.5±0.2 


-0.16±0.02 


0.02±0.03 


0.18±0.02 






-0.05 .. 


. 0.05 


2.05 


0.00 


717 


527 


2.10±0.02 


2.47±0.02 


3.50±0.02 


-0.26±0.01 


5.3±0.2 


-0.13±0.01 


0.03±0.02 


0.16±0.01 






0.05 .. 


. 0.15 


2.06 


0.09 


415 


127 


2.25±0.02 


2.51±0.03 


3.45±0.02 


-0.32±0.01 


3.7±0.1 


-0.05±0.02 


0.09±0.03 


0.14±0.02 






0.15 .. 


. 0.25 


2.06 


0.18 


130 


261 


2.38±0.03 


2.47±0.05 


3.38±0.04 


-0.29±0.02 


3.1±0.1 


-0.02±0.03 


0.13±0.03 


0.15±0.02 


2.09- 


-2.18 


-0.25 .. 


.-0.15 


2.13 


-0.18 


304 


303 


1.94±0.03 


2.48±0.04 


3.69±0.03 


-0.30±0.02 


7.4±0.4 


-0.16±0.02 


0.02±0.03 


0.18±0.02 






-0.15 .. 


.-0.05 


2.14 


-0.09 


898 


590 


1.99±0.02 


2.51±0.02 


3.84±0.02 


-0.29±0.01 


6.7±0.2 


-O.lliO.Ol 


0.10±0.02 


0.21±0.01 






-0.05 .. 


. 0.05 


2.14 


0.00 


1331 


797 


2.03±0.01 


2.56±0.02 


3.78±0.01 


-0.28±0.01 


5.9±0.2 


-0.08±0.01 


0.10±0.01 


0.18±0.01 






0.05 .. 


. 0.15 


2.14 


0.09 


838 


696 


2.11±0.01 


2.57±0.02 


3.72±0.01 


-0.26±0.01 


5.1±0.1 


-0.05±0.01 


0.12±0.01 


0.17±0.01 






0.15 .. 


. 0.25 


2.14 


0.18 


267 


443 


2.29±0.02 


2.62±0.03 


3.53±0.02 


-0.32±0.01 


3.4±0.1 


0.04±0.01 


0.17±0.02 


0.13±0.01 


2.18- 


-2.27 


-0.25 .. 


.-0.15 


2.22 


-0.18 


254 


320 


1.91±0.03 


2.43±0.04 


4.06±0.03 


-0.31±0.02 


7.7±0.4 


-0.17±0.02 


0.13±0.03 


0.30±0.02 






-0.15 .. 


.-0.05 


2.22 


-0.09 


871 


676 


1.84±0.01 


2.54±0.02 


4.09±0.02 


-0.20±0.01 


8.5±0.2 


-0.13±0.01 


0.11±0.02 


0.24±0.02 






-0.05 .. 


. 0.05 


2.23 


0.00 


1435 


960 


1.93±0.01 


2.62±0.01 


4.05±0.01 


-0.25±0.01 


7.1±0.1 


-0.06±0.01 


0.15±0.01 


0.21±0.01 






0.05 .. 


. 0.15 


2.22 


0.09 


919 


839 


2.00±0.01 


2.65±0.01 


3.96±0.01 


-0.23±0.01 


6.1±0.2 


-0.02±0.01 


0.16±0.01 


0.18±0.01 






0.15 .. 


. 0.25 


2.22 


0.18 


264 


494 


2.09±0.02 


2.66±0.02 


3.84±0.02 


-0.18±0.01 


5.0±0.2 


0.01±0.01 


0.17±0.01 


0.16±0.01 


2.27- 


-2.36 


-0.25 .. 


.-0.15 


2.31 


-0.18 


128 


258 


1.78±0.04 


2.54±0.05 


4.37±0.04 


-0.27±0.02 


9.4±0.6 


-0.14±0.03 


0.16±0.04 


0.30±0.03 






-0.15 .. 


.-0.05 


2.31 


-0.09 


569 


636 


1.71±0.01 


2.61±0.02 


4.42±0.02 


-0.15±0.01 


10.5±0.3 


-0.13±0.01 


0.15±0.02 


0.28±0.02 






-0.05 .. 


. 0.05 


2.31 


0.00 


944 


890 


1.79±0.01 


2.67±0.01 


4.40±0.01 


-0.18±0.01 


8.8±0.2 


-0.06±0.01 


0.20±0.02 


0.25±0.02 






0.05 .. 


. 0.15 


2.31 


0.09 


617 


781 


1.97±0.01 


2.72±0.02 


4.24±0.01 


-0.26±0.01 


6.4±0.2 


0.03±0.01 


0.24±0.01 


0.21±0.01 






0.15 .. 


. 0.25 


2.30 


0.18 


146 


413 


2.09±0.02 


2.74±0.03 


4.15±0.02 


-0.23±0.01 


4.9±0.2 


0.06±0.01 


0.26±0.02 


0.20±0.01 


2.36- 


-2.50 


-0.25 .. 


.-0.15 


2.39 


-0.20 


43 


160 


1.68±0.07 


2.58±0.09 


4.68±0.07 


-0.21±0.04 


11. Oil. 3 


-0.14±0.05 


0.22±0.07 


0.36±0.05 






-0.15 .. 


.-0.05 


2.39 


-0.09 


231 


152 


1.65±0.02 


2.62±0.03 


4.71±0.02 


-0.lli0.01 


11.4±0.4 


-0.13±0.02 


0.22±0.03 


0.34±0.02 






-0.05 .. 


. 0.05 


2.39 


0.00 


377 


646 


1.76±0.01 


2.77±0.02 


4.66±0.02 


-0.23±0.01 


9.1±0.3 


0.01±0.01 


0.27±0.01 


0.26±0.01 






0.05 .. 


. 0.15 


2.39 


0.08 


230 


539 


1.85±0.02 


2.80±0.02 


4.58±0.02 


-0.24±0.01 


7.7±0.2 


0.05±0.01 


0.29±0.01 


0.24±0.01 






0.15 .. 


. 0.25 


2.38 


0.18 


20 


168 


1.92±0.05 


2.85±0.07 


4.26±0.07 


-0.18±0.03 


6.7±0.7 


0.07±0.03 


0.24±0.04 


0.17±0.03 



p 

crq 

a 

a 

CO 

o 

■:. 
-■- 

r. 
-; 



< 



Note. — I e is computed in the V-band. 
^Total number of galaxies in the stacked spectrum. 
fErfcctivc median S/N of the stacked spectrum. 

|[Olll]A5007 emission line equivalent widths arc measured after subtracting off a model for the stellar population continuum. They are expressed as negative values because they represent 
emission, not absorption. 
§Stcllar population properties arc computed using emission-corrected values of H/3, assuming AH/3 — 0.7 EW([Olll]). 
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To this we add [Mg/H] and [Mg/Fe] and show how these 
quantities also vary on slices through the FP. We quan- 
tify all of these mappings by computing chi-square fits 
that express SSP age, [Fe/H], [Mg/H], and [Mg/Fe] as 
a function of 3D location in FP space. Thus the four 
classic parameters of old stellar populations can now be 
predicted purely from a galaxy's location in FP space. 
These "structural SSP" quantities can be used to esti- 
mate the star formation histories of galaxies for which 
high-quality spectra are not available. 

Having fit SSP parameters throughout FP space, we 
are then able to estimate both the mean epoch of star for- 
mation and the duration of star formation at all points in 
FP space. The additional knowledge of the [Mg/H] and 
[Mg/Fe] mappings (not considered in Paper III) proves 
to be key when interpreting the observed stellar popula- 
tion variations in terms of star formation histories. With 
these data in hand, we return to the premature trunca- 
tion model of Paper III and are now able to make the 
model semi-quantitative with estimates of how the onset 
and duration of star formation actually vary with a and 
above and below the FP at fixed a. These results are in 
good accord with the premature truncation picture. 

This paper is organized as follows. The data and bin- 
ning method used in this analysis are described briefly in 
section^ In section|31 we present the 2D family of galaxy 
star formation histories and map it onto a cross-section 
through the FP. In section HI we quantify these trends 
by providing equations for stellar population age, [Fe/H], 
[Mg/H] , and [Mg/Fe] as functions of a and A log I e . Sec- 
tion [5] goes on to intepret the observed stellar population 
trends in terms of star formation histories and presents a 
2D schematic model for the past star formation of quies- 
cent, early- type galaxies. These are used, in conjunction 
with the results of Paper III, to support the premature 
truncation model. Section |6] discusses possible physical 
mechanisms that could produce premature truncation in 
some galaxies at fixed <r, as well as the effect of dry merg- 
ing on the trends presented here. Finally, section [7] sum- 
marizes our conclusions. 

2. DATA 

The galaxy sample used in this work is the same as 
that used in Papers I, II, and III. It consists of ^16,000 
early-type galaxies f rom the SDSS Main Galaxy Sample 
(jStrauss et all I2002D in the redshift range 0.04 < z < 
0.08. The early-type galaxy selection is described in de- 
tail in Paper I. Briefly, we require galaxies to have no 
significant detectable flux in the Ha and [On] A3727 emis- 
sion lines (i.e., line fluxes are both below a 2ct detection 
threshold), to have concentrated light profiles (i.e., the 
ratio of the 90% and 50% Petrosian radii Rgo/R^o > 2.5 
in the i-band), and to have light profiles that are better 
fit by a de Vaucouleurs profile than by an exponential 
profile. The resultant sample is strongly concentrated 
on the red sequence (see Figure 1 of Paper I). Spec- 
tra for the galaxies that meet these criteria are down- 
loaded from the SDSS Data Archive ServeiQ. The var- 
ious FP paramet ers are obtained from the SDSS Data 
Release 6 (DR6: lAdelman-McCarthv et al.ll2008h Cata- 
log Archive Server (CASfl and the NYU Value-Added 

6 http://das.sdss.org/ 

7 http://cas.sdss.org/dr6/en/ 
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Figure 1. Bins used to sort and stack "similar" galaxy spectra. 
The quantity A log I e is defined as the difference between the mea- 
sured log I e of each galaxy and the best-fitting FP relation (see 
Figure 1 of Paper II) . Because A log I e is defined at fixed val- 
ues of <r and R e , it is equivalent to — A log(M ( jy 7l / Ly)- Colors 
indicate the absolute magnitude range of the galaxies to give a 
sense of sample completeness. Dark gray points indicate galaxies 
with M r < —20.1, for which our sample is volume-limited. Red 
points indicate the faint end of the sample (M r > —19.7), where 
it is missing > 50% of the volume-limited targets. Cyan points 
indicate intermediate galaxies with —20.1 < M r < —19.7. Bins 
with low values of a and A log I e are substantially incomplete and 
are populated by galaxies from the low-redshift end of our sample 
(0.04 < z < 0.05). 



Catalog dBlanton et alll2005l) fo r Data Release 4 (DR4; 
lAdelman-McCarthv et al . 1 120 061;) . as described in Paper 
II. The individual galaxy spectra typically have signal- 
to- noise S/N ~ 20 A -1 , so we sort the galaxies into bins 
based on their FP parameters, then stack the spectra 
within each bin to obtain very high S/N mean spectra 
that span the binned parameter space. 

In Paper II, we showed that the stellar population 
properties and therefore the star formation histories of 
early-type galaxies depend strongly on a and also on AI e 
residuals at fixed a and fixed R e . However, on the FP 
midplane all the stellar population properties investigated 
in Paper II showed no dependence on R e . In this work 
we simplify the galaxy structural space by marginalizing 
over i? e , such that galaxy properties are defined in the 
2D space of a and AI e only. 

The FP for our galaxy sample is shown edge-on in Fig- 
ure 1 of Paper II. Deviations from the plane are parame- 
terized as A log I e , where A log I e is defined as the differ- 
ence between the measured value of log I e for each galaxy 
and the fitted FF0, such that 

AI e =log(/ e ) - (1.16 logo-- 1.21 logiieH- 0.55). (1) 

Note that because AI e is defined at fixed a and fixed R e , 
there is a one-to-one equivalence (although a sign change) 
between Alog/ e and — Alog(Md yn /L), since M^ yn /L oc 
<j 2 R e /I e R 2 e oc a 2 /I e R e . 

To group similar galaxies, we define a regular grid 
in AI e -a space as shown in Figure Q] The right- 
hand axis reflects the equivalence between AlogJ e 
and — Alog(Mdyn/L). Within each bin, we stack all 
galaxy spectra using a c-clipping algorithm and re- 
jecting color outliers, as described in Paper I. In each 

8 Throughout this paper, we will express A log I e as "A/ e " for 
convenience, with the understanding that this is always computed 
on a logarithmic scale. 
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Figure 2. The age- [Fe/H] hyperplane for early-type galaxies. Colors and symbols show different bins in ct, as indicated. Symbol size 
indicates bins in A/ e at inte rvals of 0.1 dex, with the smallest (largest) symbols indicating the lowest (highest) AI e bins. Dashed lines are 
lines of constant ct from the Tragcr ct al. (2000a) mctallicity hyperplane. On average, high-cr galaxies have older ages and higher [Fe/H] 
than low-cr galaxies. At fixed ct, age and [Fe/H] arc anti-correlated and this relation is driven by AI e , with low (high) AI e galaxies showing 
older (younger) ages and lower (higher) [Fe/H]. This anti-correlation between age and [Fe/H] may indicate that galaxies with low AI e 
experienced shorter-duration star formation than their high A7 e counterparts at the same a (see text for details). 



stacke d spectrum, we measure the full set of Lick in - 
dices (|Worthev et al.lll994l:IWorthev fc Ottavianilll997D . 
along with estimated errors d erived from the error spec- 
tra following the formalism of iCardiel et al.l (fl998fl We 
model the mean luminosity- weighted SSP age, [Fe/H], 
[Mg/H] and [Mg/Fe], along with the statistical errors 
in each quantity, using the stellar population models 
of iSchiavonl (120071) and th e EZ_Ages code developed by 
I Graves fc Schiavonl (|2008l see this work for details of the 
modeling process). The properties of the galaxies in each 
bin, the measured Lick index values, and the derived stel- 
lar population properties are given in Table [TJ 

As always, when interpreting SSP "age" measure- 
ments, one should bear in mind that the quoted age 
values are Balmer line strength-weighted, luminosity- 
weighted mean ages. For populations with extended 
star formation histories such as galaxies, the SSP ages 
are strongly biased toward t he age of the youngest sub - 
populations of galaxies (e.g.. lTrager fc Somervillej[2Q09j) . 
Furthermore, the SSP ages presented in this paper are 
averaged over many similar galaxies. Thus, if one galaxy 
bin has a younger SSP age than another, this does 
not mean that all the stars in all of those galaxies are 
younger. Rather, it should be interpreted to mean that 

9 All H/3 line strengths have been corrected for weak H/3 emission 
infill, following the procedure described in Appendix A of Paper 
III. 



at least some (possibly quite small) fraction of the total 
star formation in those galaxies occurred up until more 
recently in the "younger" object than in the "older" one. 

3. THE 2D FAMILY OF STELLAR POPULATIONS IN 
EARLY-TYPE GALAXIES 

3.1. The Age-[Fe/H] Hyperplane 

The first result of this paper is shown in Figured which 
plots the age and [Fe/H] results from the stellar popula- 
tion analysis for the various bins in a and AI e . Colors 
and symbols give the a values for each bin of galaxies, as 
indicated in the figure key. Symbol sizes indicate bins in 
A/ e , with the largest and smallest symbols representing 
the highest and lowest values of A7 e , respectively. A ver- 
sion of this figure was presented in Paper III to demon- 
strate that galaxy star formation histories vary smoothly 
through FP space and to discuss M+/L modeling. Here, 
we present a simplified form, which collapses all values of 
R e into a single point in order to show the main trends 
more clearly. The well-known trends of age and [Fe/H] 
with a are recovered: age and [Fe/H] both increase with 
increasing a. 

Furthermore, galaxies within a given a bin span a 
range in age and [Fe/H], with age and [Fe/H] anti- 
correlated at fixed a such that older galaxies have lower 
[Fe/H] than their younger counterparts at the same o F 7 ! . 

10 Although the anti-correlation between age and [Fe/H] at fixed 
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Figure 3. (a) [Mg/H] as measured in the stacked spectra, shown as a function of [Fe/H], As in Figure[2] colors and symbols show different 
bins in cr, while symbol size indicates bins in AI e . At fixed cr, both [Fe/H] and [Mg/H] increase with increasing AI e . The slope of the 
relation is shallower than one-to-one (dotted line), indicating that [Fe/H] varies more at fixed a than does [Mg/H], [Mg/H] is always higher 
than [Fe/H], indicating that all galaxies have super-solar Mg/Fe ratios, with an offset that increases with a. (b) The same plot for [Mg/Fe]. 
At fixed <r, [Fe/H] and [Mg/Fe] are anticorrelated with the strongest anticorrelation in the highest-cr galaxies. Galaxies with low [Fe/H] 
span a range of [Mg/Fe], but all Fe-rich galaxies have relatively low values of [Mg/Fe]. The lowest A/ e bins for the three lowest values of 
a (outline d wi th black circles) deviate significantly from the [Mg/Fe] values expected from the rest of the data, and are discussed in detail 
in section [5,41 



This anti-correlation of age and [Fe/H] is essentially 
identica l to the "metallicity hyperplane" of iTrager et al.1 
(|2000al ) except that our iso-cr lines are slightly c onverg - 
ing and curved, whereas those of ITrager et al.l (|2000al) 
were straight and parallel (hence their term hyperplane). 
Lines of constant a from the Trager et al. Fe-hyperplane 
are overplotted as dashed lines. As expected, these run 
nearly parallel to our SDSS data for a g iven value of 
a. The re is a zeropoint offset between the ITrager et al.l 
(2000a) hyperplane and our iso-cr lines due to zeropoint 
differences between the stellar population models used in 
the two analyses. With the updated models used in this 
analysis, we do no t find galaxie s with a ges older than 12 
Gyr, whereas the ITrager et al.1 (|2000a| ) sample contains 
galaxies with ages > 18 Gyr. 

An important result of this series of papers is that the 
anti-correlation of age and [Fe/H] at fixed a is driven 
by variations in A/ e , or equivalently by variations in 
Alog^Mdyn/L). This was shown qualitatively in Paper 
II and will be quantified below. The correlation between 
galaxy star formation histories and structural properties 
is such that the oldest, most Fe-poor galaxies in a given 
a range are the galaxies with the lowest values of A/ e , 
while the youngest, most Fe-rich galaxies are those with 
the highest values of AI e . Variations in A/ e are evidently 

a is similar to the direction of the correlated errors in the stellar 
population modelling, these correlated errors cannot be responsible 
for the observed trends. Firstly, the errors from the modelling 
(shown in Figure [2} are too small to contribute substantially to 
the observed spread in the data. More importantly, if the spread 
in the data were due to random errors in the stellar population 
modelling, it would not be tightly correlated with A7 e . This is 
discussed in detail in Paper 1 (see section 5.4 and Figure 13) and 
Paper If (see section 4.3 and Figure 11). 



associated with differences in star formation history be- 
tween early- type galaxies with the same a. 

In Paper III, we demonstrated that stellar population 
M+/L effects can only account for a small fraction of 
the variation in I e through the thickness of the FP. This 
means that variations in I e must be due primarily to 
variations in the effective stellar mass surface density 
in galaxies (£*), assuming that our constant-IMF and 
model stellar mass estimates are valid. Thus the oldest, 
most Fe-poor galaxies at a given value of a are those with 
the lowest £*. 

The directions of increasing cr, I e , £*, and Md yn /L 
are indicated on Figure [2] by black and gray arrows, as 
labelled. These arrows reveal that the 2D family of stellar 
population properties can be mapped onto a cross-section 
through the FP, with a along one axis and AI e along 
the other. In other words, a galaxy's location above or 
below the FP is related to how it formed stars, and the 
metallicity hyperplane shows how star formation history 
maps onto slices through FP space. This mapping is 
explored further in M3 .31 

There have been some indications of this varia t ion in 
prev ious work. As menti o ned a bove, iForbes et al.l ([1998D 
and iTerlevich fc Forbes! (|2002| ) found that galaxy ages 
correlate with residuals from the FP, but in neither case 
do they identify or discuss the correlation of FP resid- 
uals with galaxy metallicities. Simultaneous with our 
Paper II, iGargiulo et al.l (|2009l ) reported similar results 
based on 141 early-type galaxies in the Shapley super- 
cluster, which show that stellar population properties, 
and [Mg/Fe] in particular, correlate with residuals from 
the FP. Their observational results agree with the results 
of Paper II, although they parameterize FP residuals in 
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terms of R e residuals at fixed a and fixed I e and therefore 
present the correlations differently. 

3.2. Trends in Mg Enrichment 

Like [Fe/H], [Mg/H] also shows variation with both a 
and with A/ e but the detailed behavior is different. Fig- 
ure [5^ shows [Mg/H] versus [Fe/H] for the various bins 
in a and AJ e . Colors and symbols are as in Figure [2] 
and solid lines connect bins with the same a. At fixed 
a, [Fe/H] and [Mg/H] are correlated with a slope some- 
what shallower than the one-to-one relation (shown as 
the dotted line for reference). The shallowness of the 
slope indicates that [Mg/H] changes less at fixed a than 
does [Fe/H]. 

The slope of the relation appears somewhat steeper for 
the lower three a bins than it does for the higher three 
a bins. However, this effect is almost entirely due to 
three data points: for the lower three a bins, the lowest 
AJ e points (outlined with black circles) have significantly 
lower values of [Mg/H] than would be expected from a 
simple extrapolation of the relation for the remaining 
AJ e bins at that a. We return to these outliers below. 

Galaxies at all values of a are offset vertically from the 
one-to-one relation (dotted line) such that [Mg/H] is al- 
ways higher than [Fe/H]. This indicates that the galax- 
ies in this sample have super-solar Mg enhancements. 
Futhermore, the size of the offset increases with increas- 
ing <7, indicating that the Mg enha ncement scales with a. 
This result is by now famili ar (e.g.. | Kuntschner fe Daviei 
[1991 Ij0rgensenlll999l [Trager et al.ll2000al) . 

Another way to parameterize the Mg-enrichment pro- 
cess in galaxies is to compare Mg abundances with those 
of Fe by looking at [Mg/Fe]. This abundance ratio is 
not a new independent variable — it is in fact defined as 
[Mg/H] — [Fe/H] — but it is conceptually useful in the fol- 
lowing discussion to examine Mg trends both in terms 
of [Mg/H] and [Mg/Fe]. In Figure |3b, we plot [Mg/Fe] 
as a function of [Fe/H] for the stacked spectra. Colors 
and symbols are as before. Galaxies appear to populate 
a wedge shape in this diagram, such that all galaxies 
with high values of [Mg/Fe] have relatively low values of 
[Fc/H] , while galaxies with low [Mg/Fe] span a range in 
[Fc/H]. The galaxies with the highest values of [Mg/Fe] 
all have modest values of [Fe/H], suggesting that high 
Mg-enhancements come at the cost of high Fe produc- 
tion in chemical evolution. 

Previous investigations of the "metallicity hyperplane" 
have not included [Mg/Fe] among the st ellar popula- 
tion p r operties that vary at fixed a fe.g.. iTrager et all 
l2000at iSmithet all l2008aft . while iThomas et all (|2005f ) 
quantify the intrinsic variation in [a/Fe] at fixed a (they 
find A [a/Fe] ss 0.05 dex) but do not correlate this scat- 
ter with other galaxy observables. Thus the addition of 
[Mg/Fe] to the spread of the metallicity hyperplane at 
fixed a is new t o this series of papers (see also the simul- 
taneous work of lGargiulo et al.N2009D . 

In this space, the three outliers identified in Figure 
[3^ stand out even more clearly. For the three lowest a 
bins, the lowest-A/ e galaxies (outlined with black circles) 
have unexpectedly low values of [Mg/Fe], which lie as 
much as 0.14 dex below the value one would derive from 
an extrapolation of the other data. It is worth noting 
that, while their values of [Mg/Fe] are unexpectedly low, 
their [Fe/H] values are a reasonable extrapolation of the 
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Figure 4. Mean light- weighted age versus [Mg/Fe], as measured 
in the stacked spectra. Colors and symbols are as in Figures l2H3l 
The three o utlyi ng data points from Figure [3p stand out strongly 
(see section 15.41 . Aside from these, SSP age and [Mg/Fe] are 
strongly correlated, with Mg-enhanced populations having the old- 
est ages. 

higher-A/ e data at the same a. This suggests that it is 
the Mg enrichment process, not that of Fe, that deviates 
in these galaxies. We return to this point in section f5. 41 

Figure 0] plots the mean light- weighted ages against 
[Mg/Fe] for each of the stacked spectra. Colors and sym- 
bols are as in Figures [2H3J Ignoring the outlier data 
for the moment, there is a tight relation between age 
and [Mg/Fe], such that galaxies with the oldest stellar 
populations are the most strongly Mg enhanced. Keep- 
ing in mind that the observed SSP ages are biased to- 
ward young subpopulations, this means that galaxies 
with some quantity of more recent star formation (i.e., 
those with SSP ages < 6 Gyr) all have relatively low lev- 
els of Mg enhancement ([Mg/Fe] < 0.2). Galaxies with 
substantial Mg enhancements have uniformly old popu- 
lations. 

In this projection, the three outlying data points with 
surprisingly low [Mg/Fe] stand out very strongly: they 
do not follow the tight age- [Mg/Fe] relation of the other 
galaxies. The discrepancy varies systematically with ct, 
being largest for the lowest-cr bin. All three points de- 
viate in a consistent way, suggesting that this deviation 
is real. This means that any physical explanation for 
the strong age- [Mg/Fe] correlation must also explain why 
low-cr, low-£* galaxies might deviate from it. We will re- 
turn to these galaxies in section 15.41 

3.3. Mapping the Age-[Fe/H] Hyperplane onto a 
Cross-Section Through the FP 

It is clear from Figures [2HD that the stellar popula- 
tion properties of early-type galaxies depend on both 
a and AJ e in systematic ways. The correspondence 
between the structural properties of galaxies and their 
stellar populations can be represented by mapping age, 
[Fe/H], [Mg/H], and [Mg/Fe] onto a-AI e space. This 
mapping is shown in Figure [5] As defined, <r-AI e space 
corresponds to a cross-section projected through the FP 
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Figure 5. Typical stellar population properties of early-type galaxies, mapped through cr-AI e space. This space corresponds to a cross- 
section through the FP. Gray points show the values of a and AI e for individual galaxies in the sample. The overlaid color contours show 
the results of our stellar population analysis of the stacked spectra. Age. [Fe/H], [Mg/H], and [Mg/Fe] all increase with a. At fixed cr, age 
and [Mg/Fe] decrease with increasing AI e , while [Fe/H] and [Mg/H] increase with increasing A/ e . 



(i.e., marginalizing over R e ). In each panel of Figure El 
gray points indicate the location of the sample galaxies 
in a-AI e space, analogous to Figure [TJ The data have 
been overlaid with color contours representing the typ- 
ical values of each stellar population parameter at that 
location. These color contours are computed by plotting 
the derived stellar population parameter values for each 
stacked spectrum at the point corresponding to the mean 
values of a and AI e for that galaxy bin. We then inter- 
polate between the 6x5 grid of values from the stacked 
spectra to map the stellar population parameters across 
the space. 

Shown this way, mean luminosity-weighted age (panel 
a) and [Mg/Fe] (panel c) track one another fairly closely, 
as expected from Figure SI with lines of constant age and 
lines of constant [Mg/Fe] running diagonally from lower 
left to upper right in the maps. This means that both age 
and [Mg/Fe] have their lowest values for low-cr, high-A/ e 
galaxies and their maximum values for high-tr, low-A/ e 
galaxies. 

In contrast, lines of constant [Fe/H] and lines of con- 
stant [Mg/H] run diagonally from upper left to lower 
right, so that [Fe/H] and [Mg/H] have their lowest values 
in low-<7, low-Aie galaxies, and their maximum values in 
high-er, high-A/ e galaxies. This is another illustration 
that the family of early-type galaxy stellar population 
properties is inherently two-dimensional and further that 
Fe does not behave like Mg. 

Although there is broad agreement in the behavior of 
age and [Mg/Fe] and also in the behavior of [Fe/H] and 



[Mg/H] , the trends are different in detail. The rough cor- 
respondence in the behavior of age and [Mg/Fe] breaks 
down at low a, where [Mg/Fe] becomes nearly constant 
for all galaxies. This is expected from Figure [4] which 
shows a steepening of the age- [Mg/Fe] relation such that 
[Mg/Fe] ~ 0.15 for a substantial range of ages. As for 
the differences between [Fe/H] and [Mg/H] , the contours 
have different slopes, such that [Fe/H] shows more vari- 
ation with AI e than does [Mg/H]. 



4. QUANTIFYING STELLAR POPULATION AND 
MASS-TO-LIGHT VARIATIONS IN FP SPACE 

4.1. Galaxy Properties in 2D: as Functions of a and 

AI e 

The consistent behavior of the various stellar popula- 
tion properties throughout these FP cross-sections makes 
it possible to write each property as a functi on of a and 
AI e . We use the IDL package mpfit.pro (Mar kwardtj 
2009) to fit robust non-linear least squares curves in mul- 
tiple dimensions, taking into account variable error bars. 
A good fit for age as a function of a and AI e requires 
some curvature (as suggested by Figure |5J|, which we 
include as a quadratic term in AI e . The resulting age 
relation is 

log age = 0.82 + 0.75 logcr' - 0.85 AI e - 1.9 (A/ e ) 2 
(0.01) (0.02) (0.03) (0.2), 

(2) 
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Figure 6. Comparing the functional fits to the observed data in 
age- [Fe/H] hyperspace. The data are as in Figure [2] The solid 
lines show the fits from equations [5] and [3] color-coded to match 
the corresponding cr bins. Including a quadratic (A7 e ) 2 term in 
the age fit produces curvature that closely matches that in the 
observed galaxies. Note that the lowest-A/ e galaxies in the lower 
two <t bins have ages substantially older than those predicted by 
the fits. These are the same galaxies that show unexpectedly low 
values of [Mg/Fe] in Figure |3p (see section [5.4l for details). 

where a' = er/150 km s -1 and AI e = log/ e — (1.16 log a— 
1.21 log R e + 0.55), as defined in Equation[TJ The num- 
bers in parentheses give the lex errors on the quoted 
parameters, based on the uncertainties in the measured 
stellar population parameters (see the values in Table [1] 
and error bars in Figures [HE]) • 

[Fe/H], [Mg/H], and [Mg/Fe] can all be fit well with 
simple planar relations, which give 



[Fe/H] = - 0.08 + 0.33 log a' + 0.59 A/ e 



(0.01) (0.02) 



(0.02), 



[Mg/H] = 0.12 + 0.66 logcr' + 0.30 AI e 



(0.01) (0.03) 



(0.04), 



and 



[Mg/Fe] = 0.20 + 0.28 logcr' - 0.29 A/ e 



(0.01) (0.02) 



(0.03), 



(3) 



(4) 



(5) 



respectively. For [Mg/Fe] , the three outlying data points 
(see section I3T21 were not included in the fit. 

Figure [6] provides a visual confirmation of the quality 
of these fits by comparing the age and [Fe/H] fits to the 
observed values. The data, shapes, and color-coding are 
as in Figure [2] For each a bin, we compute a curve from 
the age and [Fe/H] fits given in equations [2] and [3] The 
curve for each bin is computed with a fixed to the mean 
value for that bin (logcr = 1.95, 2.05, 2.14, 2.22, 2.31, 
and 2.39 for the six bins, respectively) and for a range of 
AI e , with -0.2 < Alog/ e < 0.2. These are overplotted 
as solid lines, with colors corresponding to the a bins 
represented. 

The fits do an excellent job of reproducing the data, 
with the exception of two points: the lowest-A/ e bins 
for the two lower a bins. In these bins, the observed 



stellar population ages are significantly older than ex- 
pected from the fits. These are the same data that show 
discrepant values of [Mg/Fe] in Figure [3] They will be 
discussed in detail in section HT41 

The quality of the fits is quantified in Figure [71 In each 
panel, the x-axis shows the values predicted by Equations 
[2] [5] based on the mean values of a and AI e for each 
galaxy bin, while the y-axis shows the stellar population 
parameters as measured in the stacked spectra. Colors 
and symbols are as in Figures EHU When the stellar 
population parameters are fit as functions of both a and 
AI e , the resulting relations are much tighter than would 
be derived from fitting just a, as can be seen by the 
substantial spread in values at fixed cr. 

In each panel, the inset histogram shows the distribu- 
tion of offsets from the planar fits expressed in units of 
the lcr measurement errors for each data point (i.e., the 
quantity offset/ error). These can be compared to a gaus- 
sian distribution centered at zero with c = 1 (overplotted 
lines). For [Mg/Fe] and [Mg/H], the fits in equations 0}- 
[5] replicate the observations within the observational er- 
rors. The fits to age and [Fe/H] are not quite as good but 
still indicate that the fitting functions given in equations 
[2] and [3] are a reasonable representation of the data. 

Overall, the good agreement between the data and the 
fits illustrates that the fitting functions give a good rep- 
resentation of the typical stellar population properties 
for early- type galaxies as functions of a and AI e . They 
provide a new tool for estimating the stellar population 
properties of early-type quiescent galaxies for which val- 
ues of cr, i? e , and L can be measured but for which high- 
quality spectra are not available. 

4.2. Galaxy Properties in 3D: as Functions of a, R e , 

and AI e 

This paper focuses on stellar population variations in 
the 2D projected cross-section through the FP. We have 
justified dropping the third FP dimension, R e , because 
Paper II showed that the stellar population properties 
of our quiescent galaxy sample are independent of R e at 
fixed a. It is however instructive to return to 3D FP 
space briefly, which we do in this section, in order to ver- 
ify (and quantify) the lack of stellar population variation 
with R e . Furthermore, the various mass-to-light ratios, 
stellar masses, and stellar mass surface densities of our 
sample galaxies do depend on R e , and we need to be in 
3D space to model these additional quantities. 

To do this, we bin the galaxies in a 3D parameter space 
defined by cr, R e , and AI e . This binning system is de- 
scribed in detail in Paper II, in which we determine the 
mean SSP age, [Fe/H], [Mg/H], and [Mg/Fe] for each bin 
of galaxies in the 3D space. In Paper III, we further mea- 
sure dynamical and stellar masses [Mdyn and M+jmf), 
as well as mass-to-light ratios and stellar mass surface 
densities using this system. 

Here, we fit each of the stellar population parameters 
(log age, [Fe/H], [Mg/H], and [Mg/Fe]) as functions of 
cr, R e , and AI e . The procedure is similar to that in the 
previous section except that here we are fitting in four 
dimensions instead of three. Again, we restrict all fits 
to be first order in each parameter with the exception of 
age, which we allow to include a quadratic term in AI e . 
These fits are summarized in the first section of Table [2] 
labelled "stellar populations." The quoted errors include 
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S 

s 



0.4 0.6 0.8 1.0 1.2 

0.75 log ct - 0.85 A/, - 1.9 (A/ e ) a + 0.82 



-0.3 -0.2 -0.1 0.0 0.1 0.2 
0.33 log a + 0.59 A/, - 0.08 




a , | i i i i | , , i i | i i , , | I 

r (d) 

V\\ . . . , i . . , . i . , . . i , 




1 1 1 1 1 1 1 1 1 1 1 1 L 



0.3 

0.2 

0.1 x 1 

0.0 s 
-0.1 
-0.2 



-0.1 0.0 0.1 0.2 0.3 0.4 
0.28 log a - 0.29 A/ e + 0.20 



-0.2 -0.1 0.0 0.1 0.2 0.3 
0.66 log a + 0.30 A/ B + 0.12 



Figure 7. Residuals from fits to the stellar population age, [Fe/H], [Mg/H], and [Mg/Fe] as functions of a and A/ e . In each panel, the 
x-axis gives the stellar population parameter value as predicted by Equations [2]{5] based on the average value of a and A/ e in each bin of 
galaxies. The values predicted by the fits are compared to the actual values measured from the stacked spectra on the y-axis. Solid lines 
show one-to-one relations. Colors and symbols represent a and AI e bins, as in Figures [2jj4] The outliers are circled in black. The insets 
show the distribution of offsets from the fits (black histogram) in units of the la measurement errors, compared to a gaussian distribution 
with /i = and a = 1 (over-plotted line). The insets demonstrate that Equations [2}{5] predict [Mg/H] and [Mg/Fe] within the observational 
errors, but that the difference between predicted and observed values for age and [Fe/H] show slightly more scatter than can be accounted 
for by the observational errors. See text for details. 



only the formal uncertainties in the fits as estimated by 
mpfit.pro. Presented in this way, the logc column of 
the table indicates how much the parameter in each row 
(e.g., log age or [Fe/H]) changes for each 1.0 dex change 
in a . It also means that each row can be converted into 
a functional form. For example, the first row of the table 
can be rewritten as an equation for log age, such that 



log age = 0.85 + 0.73 log a' + 0.02 log R' e 
-0.92A/ e -2.9(A/ e ) 2 , 



(6) 



where R' e = R e /2.5 kpc. 

The second section of Table (labelled "mass-to- 
light ratios") gives the parameters of fits for the dy- 
namical mass-to-light ratio (M^yn/L), the component of 
Mdyn/L contributed by dark matter or IMF variations 
(Mdyn/M+jMF), and the component of Md yn /L con- 
tributed by stellar population variations (M+jmf/L). 
The third section (labelled "stellar mass") also includes 
the stellar mass (M+jmf) and stellar surface mass den- 
sity within R e (E^jmf)- The fits in these two sections 
are taken from Table 2 of Paper III, except for the M* 
fits, which are newly computed here in an identical man- 
ner to the other fits. In all cases, the "IMF" subscript 
indicates a stellar mass or stellar surface mass density 
that has been computed assuming a Chabrier IMF. The 
effect of variable IMFs, if any, would be included in the 



Mdyn/M+jMF term (see Paper III for a thorough discus- 
sion of this point). 

Notice that all four stellar population parameters are 
essentially independent of R e in that they show 0.04 dex 
or less variation per dex of variation in R e . The data span 
almost a full dex in R e , but no systematic trends in stellar 
population age, [Fe/H], [Mg/H], or [Mg/Fe] are observed. 
Meanwhile, all four parameters depend strongly on both 
a and A/ e , with values similar to those from the 2D fits 
of section 14.11 The differences between the 3D values in 
Table [2] and the 2D values in equations [5HS] are in all 
c ases consisten t with the errors quoted for the fits. 

iSmith et al.l (|2009|) have suggested that aperture ef- 
fects may affect observed trends with L or R e . The SDSS 
spectra are obtained using spectral fibers with diameters 
of 3"; the fraction of the light sampled in each galaxy 
therefore depends on the galaxy's size and distance. If 
the stellar populations in each galaxy have substantial ra- 
dial gradients, these could create spurious (or mask real) 
stellar population trends with R e . Adopting their sug- 
gested metallicity gradients (see section 2.6 of that work) 
and applying them to our observed values of [Fe/H] and 
[Mg/H] , we find that the fits in Table [5] change mod- 
estly, such that the R e parameters for [Fe/H] and [Mg/H] 
become —0.06 and —0.08, respectively (as compared to 
their uncorrected values of 0.04 and 0.02). The a and 
AJ e parameter differences are within the quoted errors. 
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Table 2 

3D Fits: Stellar Pops and Mass Parameters as Functions of a, R e , and AI e 



Range Covered by Data 


Zeropoint 


log a log R e 
0.5 dex 0.9 dex 


Ale 

0.4 dex 


(A/ e ) 2 


Stellar Populations 


log age 
[Fe/H] 

[Mg/H] 
[Mg/Fe] 


0.85 ± 0.01 

-0.09 ± 0.01 

0.11 ± 0.01 

0.20 ± 0.01 


0.73 ± 0.05 0.02 ± 0.03 
0.27 ± 0.04 0.04* ± 0.02 
0.58 ± 0.04 0.02* ± 0.04 
0.27 ± 0.03 -0.005 ± 0.02 


-0.92 ± 0.07 
0.64 ± 0.04 
0.35 ± 0.05 

-0.27 ± 0.04 


-2.9 ± 0.5 

... ± ... 
... ± ... 
... ± ... 


Mass-to-Light Ratios 


log (M dyn /L) 

log (M dyn /M i , iIM F) 

log (M+jmf/L) 


0.62 ± 0.01 
0.20 ± 0.01 
0.42 ± 0.01 


0.85 ± 0.01 0.22 ± 0.01 
0.53 ± 0.02 0.22 ± 0.01 
0.31 ± 0.01 0.003 ± 0.01 


-1.04 ± 0.01 
-0.80 ± 0.02 
-0.24 ± 0.02 


... ± ... 
... ± ... 
... ± ... 


Stellar Mass 


log M+jmf 
log s *, IMF 


10.6 ± 0.01 
3.81 ± 0.01 


1.49 ± 0.02 0.80 ± 0.01 
1.49 ± 0.01 -1.20 ± 0.01 


0.81 ± 0.02 
0.80 ± 0.02 


... ± ... 
... ± ... 



Note. — a' = cr/150 km s _1 and R' e = -R e /2.5 kpc. The units of I e are Lq pc ~ 2 . All L values are 
V-band luminosities. The "IMF" subscript on values of M* and E* indicates that they are computed 
assuming a Chabrier IMF. Any variation in the IMF would be contained within the Mdyn/M+jMF 
term. See Paper III for details. 

""Correcting these values for differences in the fraction of the galaxy sampled by the SDSS 3" fiber 
gives the values —0.06 and —0.08 f or [Fe /H] and [Mg/H], respectively, based on the radial metallicity 
gradients assumed by Smith ct al. (200^). See text for details. 



These aperture-corrected values still show very weak de- 
pendence on i? e , indicating the the stellar population 
properties of quiescent galaxies are essentially indepen- 
dent of galaxy sizes. 

As expected from the lack of stellar population vari- 
ation with R e , the stellar population component of the 
mass-to-light variation, M±jmf/ L, also shows no depen- 
dence on R e . However, the dark matter (and/or IMF 
term) does depend on R e , as do M+jmf and Y,*jmf- 
Along the FP, the masses and mass distributions of 
galaxies depend on R e but their star formation histories 
do not. 

Table [2] illustrates several other interesting points. It 
shows the expected inverse relationship between AI e 
and Mdy n /L at fixed (a, R e ) in that Md yn /L changes 
by —1.0 dex for every one dex change in AI e . More- 
over, ~80% of the change in AI e is due to varia- 
tions in the IMF/dark matter component of the galaxy 
(Mdyn/M±jMF), with only ~20% due to stellar popula- 
tion variations (M+jmf/L). This is discussed in detail 
in Paper III. 

This dominance of Mdyn/M+jMF over M+jmf/L has 
a further implication: most of the change in AI e is due 
to a real change in stellar surface density, Ti+jmf- This 
is also evident from Table [2j which shows that £*jji/f 
changes by ~0.8 dex for every 1.0 dex change in AI e . 
Thus the observed AI e variations about the FP are pre- 
dominantly due to differences in the surface mass density 
of stars within R e , not to stellar population mass-to-light 
ratio effects. It is therefore not possible for galaxies to 
evolve through the FP simply by fading — their stellar 
mass densities are ab initio different. We return to this 
key point in section 15.21 

4.3. Parameterizing FP Residuals in Terms of R e 
versus I e 



Concurrently with our Paper II. lGargiulo et al. (2009) 
published a study of stellar populations in 3D FP space 



based on a sample of 141 galaxies in the Shapley Super- 
cluster. In that work, they parameterize residuals from 
the FP in terms of R e rather than I e . As in Paper II, 
they find that residuals from the FP correlate with stel- 
lar population ages and [Mg/Fe] . The choice of R e as the 
dependent parameter is based upon the fact that the FP 
was originally proposed as a tool for measuring distances 
to galaxies, and therefore the distance-dependent quan- 
tity R e was historically fit as the free parameter against 
the distance- independent quantities a and I e . 

We prefer to parameterize residuals from the FP in 
terms of I e for two reasons. The first is that compar- 
ing galaxies at fixed values of er and R e means that 
they also have the same total dynamical masses, since 
Mdyn oc a 2 R e . More importantly, along the FP mid- 
plane the stellar populations of galaxies are observed to 
be independent of R e . This is true over almost one dex 
of variation in R e , which represents a much larger range 
of variation than the limited thickness of the FP. For 
understanding stellar populations, it is therefore logical 
to treat R e as an "inert" variable in the stellar popula- 
tion modelling and to parameterize the residuals of the 
FP and their correlated stellar population variations in 
terms of I e . Of course, for distance studies using R ei it 
is appropriate and necessary to set R e as the dependent 
variable. 

5. THE STAR FORMATION HISTORIES OF EARLY-TYPE 
GALAXIES 

Having fully quantified the mapping between early- 
type galaxy structure and stellar populations, let us ex- 
amine the trends illustrated in Figures [2]-[5] in greater 
detail. The goal is to interpret the observed variations 
in stellar population properties as differences in the star 
formation histories of the galaxies. 

In the following discussion, we will differentiate be- 
tween two effects that can modify the abundance pat- 
terns of galaxies. The first is the effective yield of metals 
during star formation, due to the chemical enrichment of 
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the interstellar medium (ISM) by supernovae. Variations 
in effective yield also include differences in the relative 
effective yields of different elements, e.g., the Fe-peak el- 
ements predominantly produced by type la supernovae 
(SNe la) versus the a-elements predominantly produced 
by type II supernovae (SNe II). The effective yield de- 
pends on the IMF, such that top-heavy IMFs with larger 
fractions of massive stars produce higher effective yields 
overall, as well as higher [a/Fe] ratios. The effective yield 
also depends on gas inflow and outflow during star forma- 
tion; closed box models for star formation produce differ- 
ent effective yields than those with substantial outflows 
of enriched material or inflows of primordial-abundance 
gas. In addition, the effective yield of SNe la elements 
(e.g., Fe) that is incorporated into stars (and hence is 
measured by absorption line studies) depends on the 
time-scale for star formation, as discussed below. 

A second quantity with relevance for galaxy stellar 
abundance patterns is the baryon "conversion efficiency" 
of star formation, which measures the fraction of gas in 
a galaxy that is ultimately turned into stars. Galaxies 
with higher conversion efficiencies turn a larger fraction 
of their baryons into stars than those with lower con- 
version efficiencies. All other things being equal, this 
additional processing of baryonic material results in the 
production of more heavy elements over the star-forming 
lifetime of the galaxy and should result in higher abun- 
dances for all metals. 

5.1. a -Dependent Variations in Star Formation History 

The various stellar population trends with a have been 
known for some time and their implications discussed by 
many authors. We briefly review some of this discussion 
here before moving on in section [5721 to focus on the new 
results of this work: the variations at fixed a. 

As mentioned above, the effective yield of certain ele- 
ments may depend upon the timescale for star formation. 
Short timescales for star formation are often discussed 
in the literature as a possible mechanism for producing 
low-[Fe/H], high-[Mg/Fe] galaxies due to the time-delay 
between the production of Mg (and other a-elements) 
by SNe II and t he slower production of Fe by SNe 
la (|Tinslevl 119791 : iGreggio fe Renzinil 11983ft. Numerous 
previous authors (e.g., Worthev et alj 11993 IMatteuccil 



lr99ltlTrager ■ et al.ll2000at iThomas et afll2005| ) have sug 
gested that the duration of star formation (Atgjr) may 
be shorter for more massive galaxies, resulting in an in- 
crease in [Mg/Fe] with increasing a. Shorter duration 
star formation should result in fewer Fe-peak elements 
being incorporated into new generations of stars. This 
leads to lower leve ls of Fe relative to M g and therefore 
to higher [Mg/Fe]. IThomas et al.1 (|2005) have quantified 
this using a simple closed-box chemical evolution model 
to estimate a quantitative relationship between [a/Fe] 
and AtsF based on the chemical evolution models of 
IThomas et ail ([19991 ) . 

All other galaxy properties being equal, this model 
for effective yield variations predicts higher [Mg/Fe] and 
lower [Fe/H] for galaxies with shorter At$F- The pre- 
dicted trend an anti-correlation between [Fe/H] and 
[Mg/Fe] — is observed at fixed a in Figure U2 However, 
both [Fe/H] and [Mg/Fe] are observed to increase with 
increasing a. This suggests that another effect must also 
be at work along the u-sequence, a different mechanism 



that modifies the overall effective yield or changes the 
conversion efficiency of star formation versus a. 

One oft-cited possibility for this second mechanism in- 
volves relaxing the assumption of a closed-box model 
to include substantial quantitites of supernova feed- 
back. Galaxies with lower values of a have shal- 
lower gravitational potential wells and are therefore 
more susceptible to SN feedback, limiting the chemi- 
cal enrichme n t of the interstellar med i um by SNe (e.g. , 
Larson! Il97l lArimoto fe Yoshiil Il987l: IMatteuccil fl994t 



Pipino fe Matteuccil 12004 ) . This process lowers the ef- 



fective yield and tends to produce lower metallicities in 
lower-cr galaxies. It should therefore result in lower val- 
ues of both [Mg/H] and [Fe/H] , as seen in galaxies with 
low a. 

If both of these processes vary with a such that higher- 
er galaxies have short AtsF and are less affected by SN 
feedback, it is possible to explain the observed trends in 
[Fe/H], [Mg/H], and [Mg/Fe] with a. Changing At SF 
with a can also produce d the observed trends in age 
(e.g., IThomas et al.ll2005l ). as galaxies which shut down 
star formation earlier will have older SSP ages. This 
dual picture — overall modulation by a and variation in 
AtsF at fixed a — is used to construct a semi-quantitative 
model in the next section. 

5.2. A Toy Model for Galaxy Star Formation Histories 

The galaxy data presented in this work also show vari- 
ations in SSP age, [Fe/H], [Mg/H], and [Mg/Fe] in a sec- 
ond dimension that correlates with AI e at fixed a (e.g., 
Figure [5]). We hypothesize that, in addition to the pro- 
posed variations in AtsF o,nd SN feedback strength as a 
function of a, there are variations in AtsF in early-type 
galaxy progenitors at fixed a. 

This scenario is illustrated in Figure [8] by a set of toy 
models. This is a concrete and semi-quantitative repre- 
sentation of the model we first proposed in Paper III. In 
this schematic representation, star formation is approx- 
imated by top-hat models with constant star formation 
rates, i.e., star formation is "on" where indicated, and 
"off" everywhere else. In this simple formulation, galaxy 
star formation histories can be represented by two num- 
bers: the duration of star formation (AtsF) and the onset 
time at which star formation begins (t on ). In the figure, 
typical galaxy star formation histories are grouped ver- 
tically as a function of a, with color-coding as in Figures 
[20 The fill pattern indicates bins in AI e , as shown in 
the figure key. The three outlying data bins (see section 
13. 2[) are excluded from the figure. 

The duration of star formation is estimated for each bin 
by converting the value of [Mg/Fe] predicted by equation 
[5] into an es timate for At.qF ■ The conversion uses equa- 
tion 4 from IThomas et all (|2005[ ) ([a/Fe] w 1/5 - 1/6 
log Atsp)- We substitute [Mg/Fe] for [a/Fe] on the as- 
sumption that Mg enrichment traces the group of a el- 
ements as a whole. The Thomas et al. conversion from 
[Mg/Fe] to AtsF is based on a closed-box chemical evolu- 
tion model where star formation histories are assumed to 
have a gaussian shape with FWHM = AtsF- The simpli- 
fied "on/off" representation in Figure |8] can be thought 
of as showing the FWHM peak of star formation in a 
gaussian model. 

The resulting epoch of star formation is centered on 
the mean age predicted for each bin by equation [2] Since 
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Figure 8. Toy models for the star formation histories of early-type galaxies, based on the age and [Mg/Fe] fits given in equa tions [2] and [5j 
(see text for details). Different bins in a are color-coded as in previous figures. The three outlying data bins (see section f3 2D are excluded 
from the figure. Galaxies with low a typically begin star formation at later times and form stars over longer timescales than galaxies with 
higher <r. At fixed <r, the duration of star formation varies, with the highest-A/ e galaxies (filled bars) experiencing the most extended star 
formation and the lowest-A/ e galaxies (empty bars) experiencing the shortest duration of star formation. The range in IStgp at fixed cr 
leads to differing mean stellar ages and differing enrichment in SN la products that are in qualitative agreement with the trends shown in 
Figure 131 

AtsF has already been fixed by the value of [Mg/Fe] 
for each bin, this effectively sets the value of t on . Of 
course, these values are SSP ages rather than a true 
mass- weighted age. This means that the ages as rep- 
resented in Figure [5] are skewed toward young values and 
later t on than is probably realistic. However, the goal is 
not to produce fully realistic, self-consistent star forma- 
tion historief(3 but rather to examine trends for how star 
formation histories vary between galaxies with different 
properties. 

If we focus momentarily on "typical" galaxies that lie 
on the FP midplane (i.e., have AI e ~ 0.0), Figure [5] is 
qualit atively consistent with Figure 10 of Thomas ct al. 
( 2005) . Galaxies with higher values of ct form stars over 
shorter lengths of time and do the bulk of their star for- 
mation at higher redshift than lower- ct galaxies. They 
follow "staged onset" star formation, in which galaxies 
with lower values of ct have l ater t on — they begi n form- 
ing stars at later times (e.g.. iNoeske et aT] 120071 ). They 



11 Producing quantitative and self-consistent models would re- 
quire the development of full chemical evolution models, and any 
quantitative conclusions would depend on the detailed shapes of 
the star formation histories chosen for the models. An interme- 
diate path would be to estimate light-weighted ages based on the 
gaussia n star formation hist ories and At$F values estimated us- 
ing the Thomas ct al. (2005) prescription. However, this method 
cannot produce self-consistent results, as can be seen from the fact 
that the predicted values of Atgp for galaxies with nearly Solar 
abundance patterns are larger than the age of the universe. 



also manifest "downsizing" ()Cowie et al.lll999| ) such that 
lower-cr galaxies continue forming stars to later times 
than higher-er galaxies. 

The main modification proposed here is the additional 
variation in At$F at fixed ct, such that galaxies with 
low AI e have shorter duration star formation than those 
with high AI e . In its simplest form, this model would 
assume that all galaxies with the same ct form the same 
total M+ but over different amounts of time (the "fixed 
total M^" model). This changes the effective yield of 
Fe, and thus the abundance patterns of galaxies with the 
same ct. The galaxies with different AI e values at fixed ct 
also have slightly different values of t on such that lower- 
AI e galaxies start forming stars sooner than those with 
high-A/ e . However, this result depends critically on the 
quantification of AtsF- 

In general, we caution the reader against drawing 
strong quantitative conclusions from this simple toy 
model, for several reasons. Firstly, the quantification 
of AtsF cannot be correct in detail, as can be seen by 
the fact that solar-scale abundances of [Mg/Fe] imply 
AtsF > 15 Gyr in the prescription of iThomas et al.1 
(2005). Secondly, the SSP values for the stellar popula- 
tion parameters are skewed from the true mass- weighted 
values. This is likely to have a substantial effect on the 
stellar ages and a weak effect on the inferred meta llicities 
and abundance ratios (jTrager fc Somcrvillc 2009). Thus 
the toy model presented in Figure [5] should be treated as 
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Table 3 

Trends in the Toy Models 



Model 


S* 


[Mg/H] 


[Mg/Fc] 


[Fe/H] 


Age 


Fixed Total Af* Model 
Premature Truncation Model 


None 


None 


+ 
+ 


- 


+ 
+ 


Observations 


— 


— 


+ 


— 


+ 



Note. — Trends give the sign of the effect toward lower AI e 



an illustration of the inferred galaxy star formation his- 
tories, not a strongly quantitative model. Any quantita- 
tive comparisons with detailed galaxy evolution models 
should be made directly with the stellar population data 
and fitting functions of section [4] rather than with these 
toy models. 

Despite being only semi-quantitative, the toy models 
serve to illustrate some general points. Qualitatively, the 
star formation histories in Figure [5] can explain nearly all 
of the trends seen in Figures [2H5] and Table [2 as we now 
explain. For galaxies at fixed a, Figure [5] shows that 
those with low AI e have older ages, higher [Mg/Fe], and 
lower [Fe/H]. This trend is consistent with a scenario 
in which galaxies at fixed a all begin forming stars at 
roughly the same time, but experience varying Atgp- 
Those galaxies which have short AtsF contain only old 
stars and have high values of [Mg/Fe] and low values of 
[Fe/H] due to the fact that their star formation ceased 
before the ISM could fully enrich in Fe-peak elements 
through SN la. In contrast, galaxies with longer At$F 
form stars until later times and therefore have younger 
mean light- weighted ages. They are still forming stars 
as the ISM becomes enriched in Fe, giving their stellar 
populations lower values of [Mg/Fe] and higher values of 
[Fe/H]. 

The model star formation histories shown in Figure |S] 
imply that the thickness of the FP represents an age se- 
quence, with the most recent arrivals on the FP settling 
onto the high-A/ e layer of the FP. The youngest ages 
are seen at low a, suggesting that the most recent addi- 
tions to the FP ar e low-mass galaxies (in agreement with 
Treu ct al. 2005b) while the highest-cr galaxies all have 
ages > 9 Gyr. The "settling" onto the high-cr end of 
the FP happened at higher z, but those galaxies also fol- 
low the trend that the highest-A/ e galaxies arrived most 
recently. 

It is tempting to see these correlations between age 
and AJ e as self-consistent stellar population properties 
wherein younger galaxies have higher I e because they 
have younger stellar populations with lower M+/L. How- 
ever, as noted in section 14.21 this qualitatively plausi- 
ble reasoning does not stand up to quantitative analysis. 
The stellar population M+jmf/L variations through the 
FP are small — a by-product of the anti-correlation be- 
tween age and metallicity at fixed (a, R e ). The observed 
anti-correlation is in the same sense as the age-metallicity 
degeneracy; just as this combination tends to preserve 
color, it also tend to preserve M*jmf/L. We have shown 
in Paper III (as summarized here in Table [2]) that the 
highest-AJ e galaxies at fixed a and R e must have con- 
siderably higher surface mass densities of stars within R e 
than their low-surface brightness counterparts. This is at 
odds with the simple assumption that all galaxies with 



the same a form the same total amount of stars. 

Thus far we have described the fixed total M+ model, 
wherein AtgF varies at fixed a but the same total mass 
of stars forms. This model can explain nearly all of the 
behavior shown in Figure [5] the increase in age, [Fe/H] , 
[Mg/Fe], and [Mg/H] with increasing a, and the varia- 
tions in age, [Fe/H], and [Mg/Fe] at fixed a. But it does 
not explain the absolute variation in either stellar mass 
surface density or [Mg/H] at fixed a. The production 
and recycling of Mg from massive stars is nearly instan- 
taneous and therefore the effective yield of Mg should 
not change much with variable AtsF, as long as the total 
amount of star formation remains the same. In fact, the 
total enrichment of [Mg/H] seems to be closely tied to the 
total quantity of star formation experienced by a galaxy 
(Table [2]). As discussed in the following sections, the ob- 
served variation in [Mg/H] at fixed a provides valuable 
additional information about the galaxy star formation 
histories precisely because it is not strongly affected by 
At SF - 

5.2.1. Premature Truncation vs. Fixed Total M* Models 

There are two possible simple top-hat models that are 
consistent with Figure [5J The first is one in which all 
galaxies of the same a form the same total quantity of 
stars but spread out over different amounts of time, so 
that the galaxies have different star formation rates. We 
have called this the "fixed total M*" model. The second 
is one in which all galaxies of the same a have the same 
constant star formation rates and therefore form different 
quantities of stars over their lifetimes depending on how 
long star formation lasts. Galaxies whose star formation 
is "truncated" earlier will form fewer stars in total than 
those that continue star formation to later times. This 
model was introduced in Paper III to fit the trend in £* 
with AI e , where it was termed the "premature trunca- 
tion" model. These two models represent the extreme 
possibilities; reality may lie somewhere in between. 

These two models make different predictions for the 
conversion efficiency of the star formation process and for 
the enrichment of [Mg/H]. A fixed total M+ model is one 
in which the efficiency of converting gas into stars is the 
same for all galaxies with the same a. Thus the total M* 
model predicts that all galaxies with the same a and R e 
will have the same M+ and the same S*. In this scenario, 
galaxies should exhibit no variation in £* through the 
thickness of the FP. Because they have processed the 
same quantity of gas into stars, they should also have 
the same [Mg/H] (i.e., they should show no variation in 
[Mg/H] through the thickness of the FP). 

In contrast, the premature truncation model predicts 
lower conversion efficiencies for galaxies whose star for- 
mation is truncated early. In this scenario, the prema- 
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turely truncated galaxies (those with lower AI e ) should 
have lower M* and lower £*. Furthermore, they have 
processed less gas into stars leading to lower produc- 
tion of heavy elements and should therefore have lower 
[Mg/H]. In both of these scenarios, shorter duration star 
formation leads to higher [Mg/Fe] and lower [Fc/H] . 

The predicted chemical enrichment pattern for these 
two versions of the toy model are summarized in Table 
02 along with the observed variations. The table shows 
that the premature truncation model fits the data better 
than a fixed total M* model. In particular, the observed 
higher levels of both £* and [Mg/H] in higher A/ e galax- 
ies provide important leverage for distinguishing between 
the two normalizations of the toy model. The addition 
of the metallicity and abundance data here strengthens 
the case for the premature truncation model proposed in 
Paper III, which was based on the observed £* variations 
alone. 

5.2.2. Premature Truncation and a Simple Closed Box 
Model for Enrichment 

We have made the qualitative argument that, at fixed 
a and R e (and therefore fixed Mdyn), the observed differ- 
ences in [Mg/H], £*, and M* imply that some galaxies 
have proceeded farther in the process of converting gas 
into stars than others. This scenario can be tested in a 
quantitative way using a simple model for star formation 
and chemical enrichment: the closed box model. 

We construct such a model by making several simpli- 
fying assumptions. (1) Galaxies with the same values 
of Mdyn reside in the same mass dark matter haloes. 
(2) The baryon fraction is constant in all haloes. (3) 
All baryons start as unenriched gas, which is then pro- 
cessed through generations of star formation into heavy 
elements. (4) We use the instantaneous recycling approx- 
imation, since we are only concerned with the evolution 
of Mg enrichment, which is dominated by SNe II. In this 
model, the observed variations in [Mg/H], E+, and M* 
at fixed Mdyn can be directly attributed to the fraction 
of gas that has been converted into stars. 

A simple closed box model for chemical evolution gives 



z M S {t) = -2/Mghi 



M gas (t) 



M, 



baryon 



(7) 



where Zm s {t) is the abundance of Mg in the ISM at 
time t, j/Mg is the yield of Mg produced by SNe II, and 



M sas /M b 



is the fraction of the initial gas still re- 



maining in the ISM at time t ()Binnev fc Tremainelll987l 
equation 9.17). We construct a single-zone closed box 
model and allow stars to form and e nrich the ISM with 
Mg. T he model uses the Mg yields of Wo oslev fc Weaver] 
(1995), assuming their higher-energy "C" series explosion 
models for more massive stars and a Salpeter slope to the 
IMF above M — 8M Q . This model is shown in Figure 
[S] as the gray diamonds and can be fit using equation [7] 
with yMg — 0.0005 (solid black line). 

This model tracks the relationship between the frac- 
tion of gas consumed in star formation and the Mg- 
enrichment of the ISM under a set of simple assump- 
tions. We can compare it to the slope of the relation be- 
tween [Mg/H] and M* that is observed among the SDSS 
galaxies at fixed Mdyn- Table [2] shows that, at fixed 
M dy „ (i-e., fixed o and R e ), [Mg/H] ex 0.35 AI e , while 
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Figure 9. Comparing the observed relation between [Mg/H] and 
M^ at fixed M^ yn with expectations from a closed box chemi- 
cal enrichment model. Gray diamonds indicate the closed box 
model, which assumes instantaneous recycling of the SNe II mate- 
rial, a S alpeter IMF slope above M = 8Mq , and the Mg yields of 
I Wooslev fc Weaver (1995). They show how the Mg-enrichment of 
the interstellar medium evolves as gas is converted into stars (i.e., 
as Mgas /Mbaryon decreases) and can be represented by equation [71 
with j/Mg = 0.0005 (solid black line). The dashed line shows the 
observed slope of d[Mg/H]/d(log M+) for galaxies at fixed M^ yn 
(see text for details) . The closed box model predicts a similar slope 
for ^Mg — -^gas/Mbaryon as the observed relation for a wide range 
of residual gas fractions. The gray shaded region and histogram 
at the right show the range of Zy[ g spanned by the data. The 
observed values of [Mg/H] correspond to closed box models with 
M g as /M b ar yon ~ 0.2 (i.e., 80% of the gas has been consumed in 
star formation). 

logM* oc 0.81 AI e . Thus the [Mg/H]-logM* relation for 
a given Md yn should have slope 0.35/0.81 = 0.43. This 
slope is shown in Figure |H] as the dashed line. It closely 
matches the slope predicted from a closed-box model over 
a large range of residual gas fractions. This indicates that 
the observed [Mg/H] variation through the thickness of 
the FP is at least plausibly consistent with the enrich- 
ment expected from the observed change in M*. 

The shaded region in Figure [5] shows the range of Zy[g 
spanned by the data, with the dotted line showing the 
mean valuq_3- The distribution of Zm s is indicated by the 
histogram at the right. In the simple closed box scenario, 
the observed [Mg/H] values are reached when 10-50% of 
the baryons are still in the gas phase, i.e., when 50-90% 
of the gas has been co nverted into stars. 

Recent stud i es by IMandelbaum et ail J2006) and 
iZaritskv et al.l (|2008f) find baryon conversion fractions 
for early-type galaxies of 12-27% and ~10-40%, respec- 
tively. However, these studies probe the scale of galaxy 
groups and clusters. In these large dark matter haloes, 
the closed box model is obviously incorrect as primeval 
gas continues to accrete. At early times while the cen- 
tral galaxy is still forming stars, some of this gas may 
make its way onto the galaxy. This infall is ex pected 
to balance star formation (e.g. JDekel et a l. 2009). With 
this type of gas infall, i? rapidly approaches the yield 
IjBinnev fc Tremaind 119871 ) and higher absolute abun- 
dances are reached sooner, at lower M gas /Mb ary on- Thus 



12 The lAnders fc Grcvcssc ( 1989]) value for the solar Mg abun- 
dance is used to translate [Mg/H] into Zj^ g . 
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relaxing the assumption of a closed box should lower the 
conversion fraction necessary to produce the observed 
levels of [Mg/H] . Furthermore, gas accreted at late times 
after star formation has ceased will contribute to lower 
the baryon conversion fraction for the halo as a whole, 
as it is not used in star formation. 

Constructing a more realistic chemical evolution model 
is beyond the scope of this analysis, but it is interesting 
to verify that the observed slope of the [Mg/H]-logM* 
relation at fixed Md yn is reasonably consistent with ex- 
pectations from simple chemical evolution models. 

5.3. Alternative Scenarios 

We turn now to other possible scenarios for reproduc- 
ing the observed trends in age, element abundance ratios, 
and stellar mass surface densities through the FP. Some 
of these were also considered in Paper III. 

5.3.1. IMF Effects 

If the IMF is not constant, appropriate variations 
might be able to produce the observed variations in 
E* and M* at fixed er and R e . There are two simple 
IMF models: (1) a bottom-heavy IMF that produces 
more low-mass stars in some galaxies, leading to under- 
estimates of E* and M* relative to the true mass, and 
(2) a top-heavy IMF at early times that produces more 
compact remnants at late times, also leading to under- 
estimates of E* and M*. As discussed in Paper III, nei- 
ther of these models can explain the full set of observed 
variations. A bottom-heavy IMF leaves the relative fre- 
quencies of SNe II and SNe la unchanged and therefore 
does not produce the observed increase in [Mg/Fe] with 
decreasing E*. A top-heavy IMF at early times does 
increase [Mg/Fe] as required by the observations. How- 
ever, it also produces an increase in the total production 
of metals, leading to increases in [Mg/H] and [Fe/H], op- 
posite to the observed trend (see Tabled . Simple models 
for IMF variation cannot simultaneously reproduce the 
observed increase in [Mg/Fe] and the decrease in total 
metal content at low AI e . 

5.3.2. Differential Dissipation in Galaxy Mergers 

In simulations of galaxy mergers, the amount of dis- 
sipation changes the profile of the merger remnant by 
rearranging the dist ribution stars and dark matter with 
respe ct to R e (e.g., iDekel fc Cox 2006; Robertson et al.l 
2006) . Variations in the degree of dissipation are there- 
fore another way to produce variations in E* and M+ 
at fixed a. However, as discussed in Paper III, the ob- 
served trends are not at face value consistent with this 
picture. Highly dissipational mergers should be more 
common at early times, when galaxies were in general 
more gas-rich. This does not naturally square with the 
observation that less dissipational merger remnants (i.e., 
galaxies with lower E*) in fact have older stellar popula- 
tions at a given a (see Tabled]). It is possible to invoke 
subsequent dissipationless merging as a way to increase 
Mdyn/Mi, in systems that form early, but the observed 
trends extend down to small galaxies with er « 100 km 
s , where fewer dissipationless mergers are expected to 
occur. 

5.3.3. Differential Efficiency of Star Formation or 
Mass-Loss 



Other possible models to produce variations in E* and 
M+ at fixed a and R e include mechanisms that modify 
the efficiency of star formation or of chemical enrichment. 

One possibility is that the low-AJ e (low-AE*) galax- 
ies form stars for the same duration of time but overall 
less efficiently, converting lower fractions of their gas into 
stars. This successfully produces the lower observed val- 
ues of [Mg/H] in these objects. It does not, however, 
produce their higher observed values of [Mg/Fe] . 

Another possibility is that the low-E* galaxies form in 
lower-mass haloes so that more gas and metals escape 
in supernova-driven winds during star formation. Like 
the low-efficiency model above, this model also produces 
lower E* and lower [Mg/H] but does not produce varia- 
tion in [Mg/Fe]. 

A third possibility is that the low-E* galaxies experi- 
ence burstier star formation, such that many supernovae 
go off at once, producing more violent winds. This gives 
lower overall [Mg/H], as observed, but predicts the wrong 
trend for [Mg/Fe] ; strong winds during bursty star forma- 
tion should preferentially remove Mg, lowering [Mg/Fe]. 

It is clear that the combined constraints of lower 
[Mg/H] and higher [Mg/Fe] in low-E* galaxies are pow- 
erful tools for discriminating between galaxy formation 
models. None of the enhanced mass-loss or lower ef- 
ficiency models discussed in this section naturally pro- 
duces both trends simultaneously. 

5.3.4. Summary of Alternative Models 

Faced with meeting both the structural and stellar 
population constraints provided here and in Paper III, 
only the premature truncation model seems able to 
match the data without forcing. The combined strength 
of these two types of information proves to have sur- 
prising constraining power and shows the importance of 
mapping the stellar population properties into FP space. 
Putting these two different kinds of information together 
has created a tight test of galaxy formation models, 
which at this point seems to leave only one standing: 
the premature truncation model. 

5.4. The Outlying Data Points: Low-cr , Low-I e Galaxies 



In section [3~2l we identified three galaxy bins that de- 
viate significantly from what one would expect from a 
smooth extrapolation of the other data. These are the 
lowest-A/ e bins (-0.25 < AI e < -0.15) for each of the 
three lowest-cr bins (1.86 < logcr < 2.00, 2.00 < log a < 
2.09, and 2.09 < logcr < 2.18), which are outlined by 
black circles in Figures |31 El El and E] We consider here 
several possible explanations for these galaxies. 

These galaxies stand out most clearly in Figure |H 
where their [Mg/Fe] values appear much lower than ex- 
pected, particularly for the lowest a bin. Their age and 
[Fe/H] values appear relatively consistent with the other 
data in Figure O but become more discrepant when cur- 
vature is allowed in the age fit of equation [21 as can 
be seen in Figure [6] To summarize, these outliers have 
ages that are too old, [Mg/H] abundances that are too 
low, and [Mg/Fe] values that are way too low compared 
to what one would expect from an extrapolation of the 
other data. The low values of both [Mg/H] and [Mg/Fe] 
suggest that these galaxies lost a larger fraction of their 
SNe Il-enriched material than did their counterparts at 
the same a with higher AI e . 
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One possibility is that these galaxies are substantially 
affected by sample selection effects. They populate the 
part of galaxy parameter space that suffers from seri- 
ous incompleteness (Figure [lj. However, is is not clear 
how this would produce the observed deviations in stellar 
population properties. The galaxy bins contain enough 
galaxies to produce adequate S/N in the stacked spectra 
(Table [T]). They are mostly populated by galaxies from 
the lower-redshift (0.04 < z < 0.05) portion of our sam- 
ple. However, the quiescent, passive galaxy population 
we consider here should not have changed significantly 
in the ~ 500 Myr between z = 0.08 and z = 0.04 (par- 
ticularly those at low AJ e , which have SSP ages > 7 
Gyr), making it difficult to explain the observed devia- 
tions through the effects of evolution within the sample. 
The higher-redshift galaxies in these bins may be erro- 
neously assigned to the bins due to measurement errors 
in a or in AI e . These galaxies should contribute spectra 
that are representative of the galaxy properties in other, 
higher-cr, higher-AJ e bins. If the outliers are caused by 
erroneous bin assignment, they should represent linear 
combinations of galaxy properties of other bins. This is 
not the case; these bins have the lowest observed values 
of [Mg/H] and combined [Mg/Fe] and age values that dif- 
fer substantially from all other galaxy bins (see Figures 
ElandH]). 

A second possibility is that, if the stellar population 
properties of satellite galaxies differ substantially from 
those of similar galaxies that are central to their haloes, 
the observed deviations in these bins could be due to a 
higher contributi on from satellite g alaxies. We use the 
group catalogs of I Yang et al.l (|2007T ) to identify galaxies 
which have the highest M* in their group and identify 
these as central galaxies, while other galaxies are clas- 
sified as satellites. The satellite fraction, which is less 
than one-third in most of the galaxy bins, increases to a 
maximum of two-thirds in the lowest-er, lowest-A/ e bin. 
However, a preliminary analysis of the stellar popula- 
tion properties of central galaxies as compared to satel- 
lite galaxies with the same a and AI e suggests that the 
two types o f galaxies have similar stellar populations (see 
also lCooper et al. 2009). Furthermore, the satellite frac- 
tion of a bin depends more strongly on a than on A/ e ; 
the highest satellite fractions are found in the three bins 
with 1.86 < logo- < 2.00 and AI e < +0.05, rather than 
the three bins with 1.86 < logc < 2.18 and AI e < -0.15. 
We reserve a more detailed analysis of the stellar pop- 
ulations of satellite versus central galaxies for a future 
paper, but note here that our preliminary work on this 
topic suggests that it cannot account for the observed 
outliers. 

A third possibility is that the outliers are due to differ- 
ences in their susceptibility to feedback. These galaxies, 
with the lowest values of AI e and low values of c, are 
exactly those galaxies which we would expect to be most 
strongly affected by feedback. They have the shallowest 
potential wells (low a) and the least-concentrated mass 
profiles (low AI e and low AS*). The outlying galaxies 
behave in exactly the way predicted for galaxies that ex- 
perience substantial mass-loss driven by SNe II feedback 
(see section l5.3.3p . showing both lower than expected val- 
ues of [Mg/H] and [Mg/Fe]. The most deviant galax- 
ies are those in the lowest a bin, which would be the 
most vulnerable to the removal of Mg-enriched material. 



Moreover, the anomalies increase at lower a, futher tar- 
geting low a and feedback as the cause. 

A problem with this line of argument is that it is not 
clear why this effect would kick in suddenly at low AI e 
(there is no indication of deviation for galaxies with sim- 
ilar a and AI e > —0.15). This could be explained if, 
in addition to having the shallowest gravitational po- 
tentials, these galaxies have the burstiest star formation 
histories. The observed low values of [Mg/Fe] in these 
galaxies require the preferential loss of Mg-enriched ma- 
terial over Fe-enriched material. Multiple short bursts 
of star formation would maximize Mg ejection through 
large numbers of simultaneous SNe II, while the more 
extended production of Fe by SNe la would be less effi- 
cient at expelling enriched gas, leaving the Fe in the ISM 
to be incorporated in the next burst of star formation. 
Combining bursty star formation histories with shallow 
potential wells could produce the dramatic transition ob- 
served between the main population of early-type galax- 
ies and these outliers. 

6. DISCUSSION 

6.1. Physical Mechanisms for Truncating Star 
Formation 

The premature truncation model appears to reproduce 
the complex ensemble of galaxy data presented in this 
series of papers. Multiple physical mechanisms have been 
proposed that could provide such a truncation scenario, 
which we discuss briefly here. 



6.1.1. 



}uenchmg in Infalling Satellites 



One possible truncation mechanism for satellite galax- 
ies is the rapid quenching of star formation when a 
satellite falls into a massive halo and is stripp ed of 
gas fe g.. IGunn fc Gottl fl97l ILea fc De Young! [l976t 
iGislerl 119761) . This explanation predicts that the pre- 
maturely truncated galaxies (those with lower A/ e ) at 
a given a should be satellites in massive haloes, and 
that the lowest-AJ e galaxies are those that were ac- 
creted into massive haloes at the earliest times. There 
is some evidence that, at fixed c, galaxies in denser 
environments have slightly older ages than those in 
the field (e.g., iTrager et aLll2000bl; iThomas et all 1200 ' 
Clemens et al.l 2006; Bern ardi et al.l 120061 : ISmith et al 



2008bb iCooper et al.l 12009ft . However, it is not clear 



that these galaxies can be associated with the low-A/ e 
galaxies — they do not necessarily also show lower metal- 
licities and higher [Mg/Fe]3- There is observational ev- 
idence that this typ e of quenching m ay proceed slowly 
in massive galaxies (|Wolf et al.l I2009D and may not be 
adequately abrupt to produce the high [Mg/Fe] seen in 

13 Several authors find no dependence of metallicity or abun- 
dance pattern on environment at fixed a ((Thomas ct al] 120051 : 
IClemens et~aTl 120061: IBernardi et~aTl I2006T) wh ile the residual 
metallicity trends at fixed cr in Cooper ct al. ( 2009) arc in the oppo- 
site direction of what is observed here for low-A/ e galaxies: galax- 
ies in denser environments are older and slightl y metal-rich com- 
pared to field galaxies at the same a. The Smith ct al. (2008b) sam- 
ple of Coma galaxies show stellar population variations at fixed a 
as a function of cluster-centric radius, such that galaxies farther out 
in the cluster have younger ages, higher [Fe/H], and slightly lower 
[Mg/Fe] than their counterparts near the cluster center. These 
trends are in the same direction as the trends we observe with AI e 
and could be related. 
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high-er galaxies with low AI e . Furthermore, the major- 
ity of our sample galaxies are central to their haloes; two 
thirds of our sample galax ies are ident i fied as the bright- 
est group galaxies in the lYang et all (|2007ft catalog of 
SDSS galaxy groups. 

The role of satellite quenching as a mechanism for pre- 
mature truncation can and should be addressed in fu- 
ture work that distinguishes between satellite galaxies 
and those that are central to their dark matter haloes, 
as well as studying these trends in A7 e as a function of 
environment. 

6.1.2. AGN-driven Quenching 

Another possible truncation mechanism is very power- 
ful feedback, in which a substantial quantity of the ISM 
is heated and removed from the galaxy potential well be- 
fore it can form stars. AGN may be ab le to heat and 
remo v e substantial fractions of t he I SM (Gra nato et al.l 
20041: IScamiapieco fe Ohl [2Q0l iDi Matteo et all 120051: 
Springel et al.l 120051 : iHopkins et al.ll2008ft . possibly with 



an additional boost to the momen tum-driven winds from 
AGN-triggered s tar formation (ISilk fe Nusserl 12010 ). 
As an example, iNesvadba et al.l ( |2007l ) present opti- 
cal and radio observations of two massive z = 3.5 
galaxies with strong radio jets and calculate that the 
mass-loading of the jet material may represent sev- 
eral times 10 10 M© of ISM. Such objects may be evi- 
dence of powerful AGN feed back in action. Interest- 
ingly, Hopkins et al. (2007a) find that simulations of 
black hole growth through major-merger-driven accrec- 
tion predict a two-dimensional "Black Hole Fundamen- 
tal Plane" (BHFP), w hich agrees well with observations 
(jHopkins et al.ll2007bft . However, black hole growth in 
the simulations is regulated by the conditions needed for 
feedback to power a pressure-driven outflow, such that 
the resulting black hole mass depends on the details of 
the local potential and is surprisingly insensitive to other 
merger properties. In this paradigm, Mbh should be 
unique for a given combination of a and R e . It is there- 
fore not obvious how AGN-driven feedback would vary in 
such a way as to produce the observed premature trun- 
cation in galaxies at fixed a and R e . 

6.1.3. Supernova-driven Quenching and Differential Halo 

Assembly 

Supernova feedback provides another mechanism for 
injecting energy into the ISM. In a model where the 
SN-driven outflow rate is comparable to the star for- 
mation rate, SN feedback may be relevant f or driving 
galactic winds and outflow s at all mass scales (|Silkll2003r 
iPipino fe Matteuccil 120041 ) , although many galaxy evo- 
lution models require some further feedback from other 
sources such as AGN in order to f ully quench star for- 
mation in massiv e galaxies (e.g.. lHatton et al.1 120031 : 
iCroton et all 120061: iBower et all 12006ft . If some galaxies 
experience stronger SN feedback for a given star forma- 
tion rate, for exam ple due to difference s in mass-loading 
in the winds (e.g.. IPipino et al1l2009bl ). or if galaxies of 
a given mass form stars at different rates (e.g., due to 
different accretion histories), the galaxies with stronger 
feedback or more rapid star formation may shut down 
star formation earlier, resulting in low AI e galaxies with 
older ages, higher [Mg/Fe], and lower total metal enrich- 
ment. 



Another possibility is that the low AI e galaxies did not 
form the majority of their stars in situ, but instead as- 
sembled hierarchically at early times from many smaller 
galaxies. These small galaxies might be more suscepti- 
ble to supernov a feedback because t hey have sha llower 
potential wells jWhite fc Reesl[l97l iDekel fc Silld ll986t 
iWhite fc Frenkl 119911: iBenson et al.ll2003f ). This process 
could shut down star formation at earlier times and lead 
to the observed older ages and higher [Mg/Fe] of their 
massive descendents, as compared to galaxies at the same 
cr with higher AJ e , which may have formed more of their 
stars in situ in massive haloes. A galaxy which formed 
in many smaller pieces would also have lower metallicity, 
both because its ISM experienced fewer generations of 
stellar processing resulting in a lower effective yield, and 
because SN feedback would be more effective at removing 
metal-enriched gas from the galaxy. 

Differences in halo assembly histories are well-modeled 
by cosmological simulations. Thus the new generation of 
semi-analytic models, which include predictions for a and 
R e based on merger simulations, may prove important 
tools for testing assembly-dependent truncation mecha- 
nism. The successful modeling of these trends probably 
requires an improved treatment of satellite gal axies in the 
new generation of semi-analytic models (e.g.. iGuo et al.l 
12010ft over previous models, which often fail to cor- 
rectly model the stellar pop ulations of satellite ga laxies 
([Font et al.H2008l:IKimm et al.l2009t IPipino et al. 2009a). 

6.1.4. Quenching in Massive Dark Matter Haloes 

Yet anot her mechanism is quenching by massive 
haloes (e.g.. iBirnboim fc Dekell [20031: iKeres et all [20051: 
ICattaneo et al.l 120081 ). when haloes pass over a critical 
mass threshold (M cr it) and shock-heated gas accreting 
onto them can no longer cool efficiently. This process ap- 
pears to be able to quenc h star formation in ha loes with 
M ha i ~ 1O 12 -1O 13 M (IBirnboim et all [2007ft . When 
supplemented with gravitational heating from continued 
cosmological gas accretion, this effect may be able to 
truncate star formation in more massive haloes as well 
(|Dekel fc Birnboim[|2008f ) . In this scenario, the observed 
variations through the FP could be due to stochastic 
variations in the halo mass assembly history. If galax- 
ies with the same a and R e today exist in haloes that 
passed over M cr u at different times, or if mass was ac- 
creted differently onto these haloes, affecting the cooling 
rate, these galaxies could have different star formation 
histories. The low-AJ e galaxies, which were truncated 
earlier, would be those that crossed this critical halo mass 
threshold at earlier times. This scenario should also be 
quantitatively testable with semi-analytic models. 

6.1.5. Summary of Truncation Mechanisms 

These various scenarios for premature truncation make 
predictions about how a galaxy's status as a central ver- 
sus a satellite, the local environment, central black hole 
mass, dark matter halo mass, and mass assembly history 
should vary at fixed a for galaxies with different star for- 
mation histories and in different parts of FP space. A 
combination of further observational programs and im- 
proved semi-analytic models for galaxy formation (that 
can track values of a and R e for galaxies in a cosmologi- 
cal context) should be able to distinguish between these 
different truncation mechanisms. The analysis presented 
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here makes it possible to identify galaxies that have expe- 
rienced premature truncation, providing a powerful tool 
for exploring the quenching of star formation in galaxies. 

6.2. The Effects of Dry Merging 

The build-up of L ~ L* red ga laxies since z ~ 1.5 
(|Bell et alj|2004t iFaber et aLll2007l) appears to require a 
contribution from dissipationless "dry" mergers between 
galaxies that have already qu enched star forma tion and 
occupy the red sequence (e.g- lFaber et alJl2007t l. In this 
section, we explore the effects of dry merging on the lo- 
cation of galaxies in FP space and the observed stellar 
population trends. We consider the two extremes of dry 
merging: equal-mass major mergers between identical 
progenitor galaxies and the build-up of mass through a 
series of minor mergers accreting onto a massive central 
galaxy. These cases bracket the range of merger mass 
ratios. 

The scenario of mass-doubling through equal-mass 
dry mergers has been explored by iBovlan-Kolchin et al.l 
(|2005l 120061 ) using N-body simulations. The models use 
empirically-motivated initial conditions for the galaxies, 
combined with iNavarro et al.l ([1997) dark matter haloes 
which are allowed to evolve acti yely through the c ourse 
of the simulation. More recently, Na ab et all (2009) have 
modeled the effects of mass-doubling through a large 
number of minor mergers in a cosmologically-motivated 
hydrodynamical simulation of a galaxy which has no 
mergers with mas s ratios greater than 1:8 after z = 3. 
iNaab et al . (2009) also present a simple analytic argu- 
ment for the way in which a, R e and stellar density 
should scale in mergers, based on the virial theorem. We 
use these studies here to consider the effects of dry merg- 
ing on the observed stellar population trends throughout 
F P space. 

Boylan-Kolchin et al. (200G) find that equal mass dry 
mergers between galaxies that lie on the FP result in 
merger remnants that also lie on the FP. The merger 
remnants have velocity dispersions (<7/) similar to those 
of the progenitor galaxies (<Jj), with 07 w 1.2t7i for the 
most probable orbital parameters and at ~ I.Ocj for 
radial orbits. The remnant effective radii change more, 
with Rf rs I.6R4 for the most probably orbital paramters 
and Rf m 2.5Ri for radial orbits. This is i n reasonable 
agreem ent with the simple analytic model of Naab et al. 
(2009), who argue that equal- mass dry mergers should 
produce Of = <ji and Rf = 2Ri. 

If we assume that the progenitors are identical, lie on 
the FP, and have the stellar populations predicted for 
their location in FP space, then the resulting merger rem- 
nant will have the same stellar populations as the pro- 
genitors, roughly the same <r, a larger R e , and will lie on 
the FP midplane (i.e., it will have AI e = 0). This is en- 
tirely consistent with the observations in Paper II, which 
show that galaxies with the same a but different R e on 
the FP midplane have the same stel lar populations. The 
virial scalings of INaab et al.l ([2009D also predict almost 
no motion off the FP midplane, with (AI e )f = +0.06 
in the merger remnant, assuming the progenitors had 
(Aie)-t = 0.0. Thus equal-mass dry merging will pre- 
serve the stellar population trends identified in this series 
of papers and will not tend to move galaxies off the FP 
midplane. 

Applying the virial scaling arguments of INaab et all 



(2009) to mass doubling through minor mergers implies 
o~f = Q.5ai and Rf = 4Ri. Assuming the progenitors 
lie on the FP with (AI e )i — 0.0, the resulting remnant 
after many minor mergers should have (AI e )f = —0.12. 
Thus minor mergers will slowly move galaxies off the FP 
toward lower surface brightnesses, while also lowering a 
and dramatically increasing R e . 

If the progenitors have the stellar populations inferred 
from their location in FP space, the accretion of stars 
from the lower-er galaxies will tend to lower the SSP 
age of the galaxy, while also lowering [Mg/Fe], [Fe/H], 
an d [Mg/H ], However, the hydrodynamical simulations 
of INaab et al.l (2009) showed that the majority of stars 
accreted from small satellites wind up at large radius. 
At z = 0, accreted stars account for ~ 1/7 of the stars 
inside 1 kpc, beyond which the relative contribution of 
accreted stars slowly increases. This means that stel- 
lar populations sampled by the SDSS spectral fibers are 
dominated by the in situ component of stars, not the ac- 
creted stars. Thus the stellar population of the merger 
remnant as observed through the SDSS fiber will not 
change dramatically from that of the massive progenitor 
galaxy. 

This means that the net effect of a series of minor merg- 
ers will be to move the galaxy to lower a and higher 
R e , while hardly changing the central stellar population. 
From the virial scaling argument above, these galaxies 
will lie slightly off the FP at lower AI e . Because the 
main progenitors started out with higher er, they should 
have older ages and higher [Mg/Fe] than other galaxies 
with the same a / on the FP midplane. This is in the same 
direction as the observed trends at fixed a to lower AI e . 
However, they should also have higher [Fe/H] and higher 
[Mg/H], which is in the opposite direction from the ob- 
served trends. If we relax the aperture assumption and 
allow the accreted low-mass galaxies to contribute sub- 
stantially to the light sampled by the SDSS fiber, they 
will tend to contribute stars that are younger with lower 
[Mg/Fe], lower [Fe/H], and lower [Mg/H]. This improves 
the match to the observed [Fe/H] and [Mg/H] values, but 
washes out the trends with age and [Mg/Fe]. 

All in all, it appears that modest quantities of dry 
merging will leave the observed stellar population trends 
through FP space essentially unchanged. Very large 
numbers of minor mergers may tend to smear out some 
of the observed trends. 

7. CONCLUSIONS 

Papers I and II presented a detailed analysis of the stel- 
lar population properties of quiescent early- type galaxies. 
These population properties were mapped through two 
familiar early-type galaxy scaling relations: the color- 
magnitude relation (Paper I) and the Fundamental Plane 
(Paper II). Paper III supplemented that work with a de- 
tailed study of galaxy mass-to-light ratios, showing that 
stellar population mass-to-light ratio effects do not con- 
tribute the majority of the observed variation in AI e . 
These works illustrated that the stellar population prop- 
erties (and by inference the star formation histories) of 
these galaxies form a two-parameter family. Further- 
more, this two-parameter family maps onto the large- 
scale structural properties of the galaxies. 

In this work, we have presented what we feel is the 
clearest and most intuitive representation of the rele- 
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vant structural properties for early-type galaxies: a and 
AI e . We have shown how the various stellar popula- 
tion parameters map onto this structural space, which 
corresponds to a cross-section projected through the FP. 
We have also explored in detail the co- variation of mean 
light-weighted age, [Fe/H], [Mg/H], and [Mg/Fe] among 
early-type galaxies. By combining these stellar popula- 
tion results with the mass information from Paper III, we 
have presented a plausible scenario for early-type galaxy 
star formation histories that is consistent with the multi- 
dimensional properties of the data. Our main conclusions 
can be summarized as follows: 

1. The stellar populations of quiescent early- type 
galaxies form a 2D fa mily, similar t o the age- 
metallicity hyperplane of iTrager et a l. (2000a]). In 
the hrst dimension of the parameter space, mean 
light- weighted age, [Fe/H], [Mg/H], and [Mg/Fe] all 
increase together. In the second dimension of the 
parameter space, age and [Mg/Fe] increase while 
[Fe/H] and [Mg/H] decrease. 

2. This 2D family of stellar population properties can 
be mapped onto a cross-section projected through 
the FP, parameterized in this work by a and AI e . 
The first dimension of the parameter space maps 
onto a, such that age, [Fe/H], [Mg/H], and [Mg/Fe] 
all increase with increasing a. The second dimen- 
sion maps onto AI e at fixed a, such that as AI e 
increases, [Fe/H] and [Mg/H] increase while age 
and [Mg/Fe] decrease. 

3. Given the systematic behavior described in conclu- 
sion 2, it is possible to estimate the stellar popu- 
lation properties of a quiescent early-type galaxy 
based on its location in FP space. Equations [5]-[S] 
quantify this mapping. They are shown to repro- 
duce the binned observed stellar population prop- 
erties to roughly within the observational errors. 

4. Paper III showed that variations in AI e are primar- 
ily due to variations in stellar mass surface density, 
AY,+jmf- This means that the observed stellar 
population variations through the thickness of the 
FP correlate with variations in AS^/mf, or al- 
ternatively with variations in AMjyn/M+jMF- At 
fixed a, galaxies with younger ages, lower [Mg/Fe], 
and greater total metal enrichment (i.e., higher 
[Fe/H] and [Mg/H]) also have larger Y,+jmf and 
formed more stars. 

5. The observed stellar population variations can be 
interpreted in terms of star formation histories 
where the start time and duration of star formation 
both depend on a (e.g., a "staged star formation" 
model), but in which there are also systematic vari- 
ations in the duration of star formation (AtsF) at 
fixed a. 

6. At fixed a, galaxies with longer duration star 
formation have higher H+jmf and higher total 
M+jmf at fixed R e . They also are more en- 
hanced in both Fe and Mg, with lower [Mg/Fe] than 
their counterparts that experienced shorter dura- 
tion star formation. This strongly suggests that 



they have experienced higher "conversion efficien- 
cies" of turning their gas into stars. 

7. The observed pattern, in which longer duration star 
formation corresponds to the production of more 
stars, while shorter duration star formation pro- 
duces fewer stars, leads us to propose a premature 
truncation model for quiescent galaxy star forma- 
tion histories. In this model, star formation pro- 
ceeds similarly for all galaxies of the same a but 
is truncated at different times. Thus some galax- 
ies are shut down while their formation process is 
less complete, leaving them with fewer total stars, 
lower Y>+jmf, and lower metal abundances. 
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